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Discrete surfaces of constant mean curvature 

(Fran Burstall, Udo Hcrtrich-Jcromin, Wayne Rossman & Susana Santos; 17 April 2008) 



Abstract. We propose a unified definition for discrete analogues of constant mean curvature sur- 
faces in spaces of constant curvature as a special case of discrete special isothermic nets. Backlund 
transformations and Lawson's correspondence are discussed. It is shown that the definition gen- 
eralizes previous definitions and a construction for discrete cmc surfaces of revolution in space 
forms is provided. 
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1. Introduction 

66" , ^ 

f^ ■ Discrete surfaces of constant mean curvature (discrete "cmc surfaces" ) have been studied in re- 

^^ \ cent years from a variety of different points of view. Two essentially antithetic approaches, one 

Cn ■ from variational principles and the other from integrable systems, lead to substantially different 

5-H I definitions: normally, for example, discrete soap films and bubbles, i.e., "discrete variationally 

cmc surfaces" , are triangulated whereas the definition of "integrable discrete cmc surfaces" makes 
use of special coordinates and, therefore, leads to "discrete cmc nets'"'' , i.e., quadrilateral surfaces. 
f~^ ■ Even in cases where it is sensible to compare the two approaches, such as that of a discrete 

catenoid in Euclidean space, it turns out that different notions are obtained: each approach leads 
to a different class of discrete surfaces that can be viewed as analogues of the smooth catenoid. 

\^ • The present paper is concerned with the integrable systems approach to discrete cmc surfaces. 

A key feature of this approach is its compatibility with the transformation theory of the (smooth) 
surface class under consideration: for a given class of surface, not only is a similar transformation 
theory sought for the discrete case but a discrete surface in the class should be created by repeated 
(Backlund-Darboux-)transformations of the smooth class; or, otherwise said, every 2-dimensional 
subnet of a multidimensional net created by repeated transformation of a discrete surface in the 
class should itself be a discrete net of the class. This is what has recently been coined "multidi- 
mensional consistency" , see the very clear and essential description of integrable discretization in 



< 






f^ I A key idea in the definition of integrable discrete cmc surfaces has been to consider them as special 

[~^ . discrete isothermic nets: that is, to discretize a (conformal) curvature line net on a smooth cmc sur- 

04 ' face — recall that smooth surfaces of constant mean curvature (in any space form) are isothermic, 

^^ . i.e., allow a parametrization by conformal curvature line parameters. This has been the pioneering 



f^ ' idea in [3] , where the authors introduced the notion of discrete minimal surfaces in Euclidean space 

QP ■ alongside the notion of discrete isothermic surfaces and their Christoffel transformation . Subse- 

^""^ ' quently, the notion of discrete surfaces of constant mean curvature in Euclidean space has been 

introduced alongside a notion of a Darboux transformation for discrete isothermic nets in [16], see 
also |4|, and the notion of discrete horospherical surfaces in hyperbolic space — as an analogue of 
smooth cmc 1 surfaces — h as b een introduced alongside a notion of a Calapso transformation for 



^ ^ 

Cu ^ discrete isothermic nets in [17] . In all three cases the constant mean curvature surfaces can be 

characterized as isothermic surfaces with a special behaviour of their transformations, as we will 
discuss below. 

Note that for all three classes of surfaces H^ + k> 0, where H is the mean curvature of the surface 
and K is the ambient curvature. It is straightforward to use the Calapso transformation for discrete 
isothermic nets to extend this family of definitions to discrete analogues of any constant mean 
curvature surfaces with 

H^ + K>0. 

Here, the key observation is that, for smooth constant mean curvature surfaces in space forms, 
the Calapso transformation becomes a conformal variant of the Lawson correspondence and that 
Bianchi permutability can then be used to carry over the characterization of cmc surfaces in 
Euclidean space to other space forms. However, these ideas turn out to be useless for surfaces 



1) Note the parallel with Christoffel's original paper [12]| , where his transformation is introduced — motivated by an 
observation about minimal surfaces. 
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ith 



H' 



K< 



as, for example, for minimal surfaces in hyperbolic space. 

For discrete isothermic surfaces in Euclidean space, a construction of a mean curvature function 
or rather a "mean curvature sphere congruence" was given in [4] and shown to be constant for 



discrete minimal or constant mean curvature surfaces . Note that this mean curvature function 



is defined at the vertices of a discrete isothermic net. Very recently, new ideas from [22] have led 
to substantial progress in this direction: a new definition of discrete cmc surfaces in Euclidean 
space relies on the requirement that a mean curvature function — defined via Steiner's formula 
on the faces of a discrete (isothermic) net — be constant, see also 
to the one via isothermic transformations, see [91. 



This definition is equivalent 




Fig. 1. Discrete minimal net in H'' 



Our mission in the present paper will be to add another definition of discrete cmc surfaces to the 
list. However, the aim is not just to promote mathematical pluralism: our definition provides a 
uniform definition of discrete cmc nets in all space forms alike — in particular, we also capture 
the previously inaccessible case of 

ij2 + K < 0. 



In f act, we define the much wid er class ' of "discrete special isothermic surfaces" based on [11 



and |[21, Def. 2.18] 



Def. 3.12| ; these come equipped with a "type number" N & N — discrete 
cmc nets in space forms will be the A^ = 1 case. Hence our definition does not only provide 
a generalization in allowing any ambient space form and value of the mean curvature, but also 
in discussing a wider class of discrete isothermic nets — and we expect it to inaugurate a new 
direction of research in the field. 

We shall start our investigation with a short discussion of discrete isothermic surfaces and their 
transformations — not only to remind the reader of some facts and to fix notations but also 
to introduce our perspective on discrete isothermic nets via loops of flat connections, which will 
be central to all that follows, see Lemma 2.5. This will set the scene for the central section 



of this paper: we shall investigate the properties of polynomial loops of parallel sections, called 
"polynomial conserved quantities" , and relations to the geometry of the underlying isothermic 
net. Excluding some degenerate cases we will arrive at the notion of "discrete special isothermic 



2) In the minimal ease the reverse is in fact also true as shown in |3|. 

3) In contrast to the generic terminology of "special discrete isothermic nets" used earlier our "discrete special isothermic 
surface" will be a technical term. 
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nets of type N^^ in a natural way. It turns out that the Darboux transformation for discrete 
isotherniic nets behaves nicely on these special isothermic nets, which gives rise to a "Backlund 
transformation" for special isothermic nets; in particular, we will prove a Bianchi permutability 
theorem that establishes "SD-consistency" for special isothermic nets. Hence our discrete special 
isothermic nets satisfy the two fundamental discretization principles of the "discrete Erlanger 



programme" of [6] 



Transformation group principle — this is built into our construction as we are working in 
conformal geometry which is the na tural symmetry group for special isothermic (smooth) 



surfaces and (discrete) nets alike, see |[21, Sect. 2.2.3 



Consistency principle — which is established by our Bianchi permutability theorem for the 



Backlund transformation, see Thm. 4.7 



Certain (very) special Backlund transforms of a special isothermic net, its "complementary nets", 
will provide the basis for establishing the relation of our approach with the previous ap proac hes 
to d iscre te cmc surfaces via their transformations as discrete isothermic surfaces in |3]], |[16]| , |[4| 



and [17], as discussed above. Moreover, we obtain a characterization for discrete cmc surfaces in 



space forms, i.e., special isothermic nets of type 1, with 

via complementary nets — as one may have expected; and the lack of their existence when 

provides one possible explanation why the aforementioned approach to define discrete cmc nets in 
space forms was doomed to failure in this case. In the same context we also obtain characterizations 
of type 2 special isothermic nets, which discretize the classical "special isothermic surfaces" of 
Darboux [[is] and Bianchi [2], see also |[14, §§84-86| . 



In the final section we arrive at the main subject of this paper: we give a definition of discrete cmc 
surfaces in space forms as special isothermic nets of type 1. Clearly, the rich theory developed in 
the more general case of special isothermic nets descends to a similarly rich theory for discrete cmc 
nets in space forms — in particular, we have the discrete analogues of the Lawson correspondence 
and the Backlund transformation and our discretization satisfies the two discretization principles 
above for an "integrable discretization". Note that we consider Mobius geometry as the natural 
ambient geometry for constant mean curvature surfaces in space forms as these arise in "Lawson 
families" of cmc surfaces with different ambient curvatures: the confinement to a space form 
subgeometry appears as a symmetry breaking phenomenon initiated by part of the geometric 
data attached to a special isothermic net and, in particular, a discrete cmc net. 

Despite the obvious merits of our definition we make a great effort to convince the reader of its 
value by providing detailed analysis of how the aforementioned previous approaches tie in with our 
definition. However, we do not conceal its problems: we provide an example of a single discrete 
isothermic net which is a cmc net in a whole family of different space forms ; even though this 
seems to be a rather singular example we can, at the moment, only speculate about how to 
obviate this anomaly. On the positive side, we provide a method to construct discrete cmc nets 
of revolution for any prescribed mean curvature H and ambient curvature k. In particular, we 
show how to explicitely construct discrete analogues of smooth constant mean curvature surfaces 



that were previously unavailable: for example, we construct (see Fig. 1 ) the discrete analogue of 
a "hyperbolic catenoid" , that is, a minimal surface of revolution in hyperbolic space, see |l]. 
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U. PiNKALL, W. SCHIEF and Y. SuRiS for many interesting and helpful discussions about the 
subject. 
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The surface graphics in Figures [^ & gwere produced using Mathematica. 



4) This net is not spherical but it "looks" close to a "wrinkled" sphere — note that, in a smooth world, spheres are the 
only surfaces that have constant mean curvature in different space forms. 
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2. Discrete isothermic nets 

We consider discrete nets f : Z^ D M ^ S^ in the (conformal) 3-sphere, where 

M — {(to, n) G Z'^ \ TO-i < TO. < TO.2, Til < n < TI2} 

is a rectangular grid : 

Def. 2.1. Such a net will be called a discrete isothermic net if there is a (real) function a on the 
edges of M, that is, a map (ij) i-^ Uij e R with Uji = aij for all edges (ij), so that 

(i) a has equal values on opposite edges of elementary quadrilaterals 

(ijkl) = ((to,, n)(m + 1, n)(m + l,n + 1)(to,, n + 1)), 

j'-e-, a(„^„)(™+i_„) = a(,„_„+i)^(„+i_„+i) and correspondingly for "vertical" edges; 
(ii) the cross ratios ' qijki — [fi] fj',fk', fi] on faces factorize as 

aij 

lijkl — — 

into two functions of one variable. 



Thus we employ the "wide definition" of discrete isothermic nets |[4|, see also [18, §5.7.2] , which dis 



cretizes isothermic nets parametrized by curvature line coordinates (not necessarily conformally) : 
as all cross ratios are real, the four vertices of any face of the net are concircular, so that a discrete 
isothermic net qualifies as a discrete curvature line net or discrete principal net. Note that the 
smallest domain of a discrete net where "discrete isothermic" imposes a condition is a 3 x 3-grid, 
TO-2 — TO,i = 712 — ni — 2; there the definition can be reformulated as a cross ratio 1 condition on 
four cross ratios: 

9(m,,l-l)(m+l,„-l)(,T.+ l,,^)(n^,,^) 9(m - 1 . ,1) (,ti ,,^) (,n,,i+ 1) (,n- 1 ,,^ + 1) _ ^. 
^(m.riXm + l, „)(,„+!, „ + l)(m,„ + l) g(™ - 1,„ - 1) (m.n- 1) (m,,i) (m - 1 ,„) ' 

a cross ratio function, satisfying this condition on all 3 x 3-grids in M , determines the function a 
uniquely up to a non-zero factor. 

As a mild regularity assumption, discretizing the notion of an immersed surface parametrized 
by curvature lines, we will usually add the requirement that any three of the four vertices of a 
face uniquely determine the circle of the four vertices, i.e., that any three vertices are in "general 
position" . 

Throughout the paper we will use the following notations: if g is a map defined on the vertices of 
a rectangular grid M, then we let 

d9ij-=93-9i and gij := \{g,+ gj); (2.1) 

note that (ij) i— > g^j defines a function on the edges of M whereas (ij) h^ dgij defines a 1-form, 
that is, dgij + dgji = 0. With these notations a Leibniz rule holds: 

d{g ■ h)ij = gij ■ dhi-j + dgij ■ h^j, (2.2) 

where "•" denotes any product on the target space of g and h. 



5) It should be straightforward to generalize our results to discrete nets defined on quad-graphs, making it possible to 



consider discrete isothermic nets with "umbilics" , cf . 

^•3 



19 



6) Note tlia t tlie cross ratio of four points in S is (up to complex conjugation) a conformal invariant. For a detailed 



discussion sec 



18, Sects. 4.9, 6.5 and §7.5.14 
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2.1. The projective approach and Moutard Hfts 

As we are considering nets in the conformal 3-sphere it will be helpful to consider 

S^^L'^/RCEP'^, where L^ = {Y e R'^'^ \\Y\^ = 0}, 



as a quadric in projective 4-space. Recall (from [18] for example) that 2-spheres are, in this model, 
described by Minkowski 4-spaces in IR^'^ and circles by Minkowski 3-spaces or, equivalently, by 
their (spacelike) orthogonal complements, and that incidence translates into a subspace relation 
or as orthogonality, respectively. For example, four points Pn € S^, n = 1, . . . ,4, generically lie 
on a unique 2-sphere S C S^ which can be described as 

S'9s^span{Pi,...,P4}C JR4'\ 

where P„ £ L^ with JRP„ = Pn G S''^ = L'^/R — their (complex) cross ratio ' is given by 



q = [pi;p2;P3;pi\ = 2(Pi,p4)(P2.p3) ; (2.3) 

the cross ratio becomes real exactly when the P„ become linearly dependent, i.e., when they span 
a 3-dimensional Minkowski subspace and the four points Pn are concircular. In this case the (real) 
cross ratio uniquely determines the relative position of the four points on the circle: given three 
of the points, say pi, p2 and p4, and a real cross ratio q the fourth point pa = F^^ p^iPi)j where 

F^,^, iX):^X + jj^{iq- 1) {X, P')P + (i - 1) {X, P)P'}, (2.4) 



as is easily verified from (2.3); note that F^ , € 0(4, 1) descends to a Mobius transformation of 



S^ and does not depend on the choice of representatives P, P' € L'^ oi p,p' g S^. Also note that 

jRpi ^ Ku {oo} 3q^ r«_p,(p") e 53 

yields a 1-to-l parametrization of the circle through three distinct points p,p' ,p" € S^ in terms 
of the cross ratio, so that T"py{p") = p', Tlp,{p") = p" and T^p,ip") = p. 

Now suppose that / : M — > 5"^ is an isothermic net and fix a cross ratio factorizing function a. 
We wish to show that there is a lift F oi f with 

(P„P,)=a,, (2.5) 

on every edge (ij) of M. To this end we have to show that this scaling is compatible on any 
quadrilateral: thus let {ijkVj denote an elementary quadrilateral and choose a light cone lift 
Fi e L^ of fi] then we normahze hfts Pj,Pz G L* of fj and /; so that {Fi,Fj) = Uij and 
{Fi, Fi) = an. Now we choose the lift 

^fe --^ ^th (^^) = ^« + WtM^Pj - Pi) (2-6) 

of /fc and readily verify that {Fj,Fk) = aa = ajk and {Fi,Fk) = a^ = au- 

Note that this lift F oi f satisfies the discrete version |(2.6) of a Moutard equation. 



Conversely, if a light cone lift P of a discrete surface satisfies a Moutard equation, Fk—Fi \\ Fj — Fi 



on all faces, then it is isothermic, see [7, Def. 9] . Namely, taking scalar products we learn that 



Fk + F,±F,-F,\ r {F,,Fk) = (P., Fi) 

Fk~F,±Fj+Fij ^ \ (Pfc, Fi) = (P„ P,-) 

and hence 

{F,,Fk){F„Fi) = (P„ F, - Fi) (Pfc - F,,Fi) = ((P,;, P,) - {F,,Fi)f 



7) Note that dct((Pi, Pj)) < so that ^<iet[(Pi,Pj)) 6 iJR. Using the Clifford algebra of Hfi_aClifibrd algebra valued 



cross ratio can be defined whose "imaginary" part encodes the 2-sphere of the four points 

5 



18, Sect. 6.5 
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SO that (2.3) gives 



_ {F.,F,) 



[ji) Jj) Jk', Jl\ — {Fi,Fi) 



From (2.6) it is also straightforward to see that any diagonal vertex star of a discrete isothermic 
net is cospherical: if i(m .n), rrijU E { — 1, 0, 1}, denote the vertices of a 3 x 3-grid then the discrete 
Moutard equation (2.6) shows that the four diagonals Fi — Fp , m,n E {±1}, are linearly 
dependent so that 

dimspan{F,(„ ^, , F,^^_^^ , Fi^_^^^^ , F,^_^_^^ , i^^^^^^, } < 4 

and the five points lie on a 2-sphere. 

Assuming that the vertex star {Fi^^ ^^ , Fi^-^ ^^ , Fi^^ ^^ , Fi^_^ ^^ , Fi^^ _^j } is not cospheric al the conve rse 



can also be shown, leading to two characterizations of discrete isothermic nets, see [7, Sect. 3] 



Lemma 2.2. A discrete net f : Z^ D M ^ S^ in the conformal 3-sphere is isothermic iff 

(i) there is a lift F : AI ^ L^ of f satisfying a discrete Moutard equation Fk — Fi \\ Fj — Fi on 
every face {ijkl) iff 

(ii) any diagonal vertex star is cospherical. 

The sphere containing a diagonal vertex star of an isothe rmic net is referred to as the central 



sphere of the net at the center of the star, see [7, Thm. 10] 



As an example we investigate discrete surfaces of revolution: consider the discrete net 

(m, n) y-^ /(„,„) := (?7„i, Qm cos ipn, g^ sin(p„) E R^ (Z R^ U {oo} ^ S^ 

where 77, g and ip are real functions of a discrete parameter. A straightforward cross ratio compu- 
tation would reveal that 



^(m,n)(m+l,n)(m+l,n+l)(m,n+l) 



_ (ri')^.,„ + i)" + (de^.,„+i)' 



4:Qm Sm+1 sin' 



2 <"fin,n + l 



identifying the net as a discrete isothermic net. However, we shall proceed differently to show 
that / is an isothermic net and to find a cross ratio factorizing function a on the edges. 

Consider 

as the quadratic form of the Minkowski scalar product of K^'^ and let 

F^:={^,L'-^) (2.7) 

denote the Euclidean lift ' of / into the light cone L^ C JR^'^. Now observe that 

F^m,n) := ^i^f™,„) = (-l)'"{( ^l±g±g^, 1^, 0, 0, '-f^-/- ) + ( O,0,cos^„,sin^„,q )}, 

= : M„ E ffi2>i - ^nC E R^ 

that is, there is an orthogonal decomposition IR^'^ = JR^'^ ® R^ so that 

Fi„,,n) = (-1)'"(M„ + $„C) = (-l)"$„(A/„ + C), (2.8) 

where $„ are rotations of JR^ , C E S^ El R? and M takes values in the hyperbolic plane ' 

H^ = {Y E JR2^i I |rp = -1, Fo > 0} C JR2'1. 



8) The Euclidean lift into the (flat) quadric of constant curvature (see [18, Sect. 1.4 



Q = {Y eL-^\{Y,Q) = -l}, where Q := (1, 0,0,0, -1). 
9) Secretly we are using a cnnfnmial map JR'' \ {axis} -^ H^ x S^ G JR^'^ ® JR^ = JR*'"' adapted to the rotational 



symmetry of the map /, cf. [18, §1.4.16 
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In particular, Mm -L ^nC for all (m, n). 

Note that M^ — span{$„C|n e Z} defines an elliptic sphere pencil, hence (cf. [18, Sect. 1.2] ) a 
circle 

L^ n jr2,i ^ l4 p |$^^(^ I ^ ^ ^|±^ 

which is the axis of our discrete surface of revolution. At the same time, it is the infinity boundary 
of the hyperbolic 2-planc H^ of the meridian curve. 

Clearly, F satisfies the discrete Moutard equation 

F(rn+l.n+l) - F(^m.n) = (-1)™+ {Mm+1 + ^n+lC + Mm + ^nC} = -F'(m+l,n) - F(m,n+1) 

and is therefore a discrete isothermic net with cross ratio factorizing function 

„^. .^ ^p^ p.^ ^ I -1 - {Mm,Mm+i) = ^'"""'""^ftl?:""^'^^ fo^ fe') = iim,n)im + l,n)) 
'' ' "' l-l+($„C,$„+iC)=-2sin2%^ forfo-) = ((m,n)(m,n + l)) 

as soon as Mm+i ^ Mm and ^n+iC ^ <f>„C. 

2.2. Quaternions and the Calapso transformation 

The Calapso transformation, or T-transformation, of (discrete) isothermic nets will be central to 
our investigations — it was introduced in |17] (see also [18, §5.7.16]) using a quaternionic setup for 
Mobius geometry. Hence we will first briefly discuss the quaternionic approach in order to make 
contact with earlier work; however, we will provide an independent definition in the following 
section so that a reader unfamiliar with previous approaches may just skip this section. 

Thus, we consider 5"^ = ImH U {oo} C HP^ and 



JR4.i^ 



1 ^ {X e End(M^) \ X = f ^ "y) ,x£lmM,xo,Xoc E R} C End(M'^) 

\ Xq X J 

equipped with \X\'^ = —X"^ = x^ -\- xqXoo as the quadratic form of the Minkowski product . In 
particular, we obtain an isometry 

iR3^ImH9a;^X= f^ :^') eL^c JR"'!; (2.9) 



1 -X 

4 



note that, for two such "Euclidean" light cone lifts X, F G L , 

-2(X, Y) =XY + YX = -{y - xf = \y - x\^. 
In this setup the Mobius group 

Mob(3) = {(c d) ^ End{]H^)\ac + ca = bd + db = 0,ad + cbe R\{0}}/R(1 PG1{2,H) 
of S^ acts isometrically on R*'^ via 

Mob(3) X JR44 3 {A,X)^ A-X -.^ AXA-^ G K^.i. (2.10) 

Now let / : M ^ iS*"^ be a discrete conformal net, with cross ratio factorizing function a on the 
edges, and define 



where F is any light cone lift of /; assuming that f : M ^ ImH C S^ and using the lift |(2.9) 
find 



we 



(''i"-i{;)' -^^'^^ df*,^a.M^,r 



10) This is analogous to the Vahlcn matrix approach to Mobius geometry (see [18, Sect. 7.1] ) using the Chfi'ord algebra 
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is the "derivative" of the Christoffel transform /* of / in R^, see |[4, Thm. T4|| or |[18, §5.7.7] 
Then, for A e JR, 

(1 + XTij){i + \Tji) = 1 - Afly e R 

and 

(1 + ATy)(l + Ar^-fc) = (1 + AT,i)(l + Xrik) 

on every elementary quadrilateral (ijkl) so that 

(ij) ^ 1 + Ar,, , 1 + An, : {j} x 5^ ^ {i} x S\ (2.12) 

defines a flat Md b(3)-connection ' on M x S"^, as long as 

1 - Xttij 7^ ^ 1 + ATy G Mdb(3) 

for all edges (ij). Hence there is a gauge transformation 

T^ : M ^ Mdb{3), T/ = T^(l + At^), (2.13) 

which identifies the (1 + AT)-connection on M x S^ with the trivial connection. 
The T play a key role in the transformation theory of (discrete) isothermic nets; in particular, 
it turns out that every T^f defines ' a discrete isothermic net: T"^ : M — > Mdb(3) are the 
Calapso transformations of / and the discrete isothermic nets T"^/ are its Calapso transforms, see 



18, §5.7.16] 



The connections {ij) i-^ 1 + Ar^j lift to fiat 0(4, l)-connections (ij) i~^ F^ on the (discrete) vector 



bundle M x JR^-i via |(2.10)| to give 

X ^ r^^. . X := ^4^(1 + At,,) X (1 + Xr,,) = T^-p {X). (2.14) 

2.3. The vector bundle approach 



Clearly, the fiat connections F'^ on M x R'^'^ in (2.14)| can be defined without reference to the 
quaternionic approach. First we define our setup: 

Def. 2.3. A connection on a (discrete) fibre bundle F -^ M , where the base M is a rectangular 
grid as before, is a map that assigns to each directed edge (ij) in M an isomorphism 

Fy : F^ -^ F, so that F„Fj, = 1; 

it will be said to be a fiat connection if its holonomies around all elementary quadilaterals (ijkl) 
are trivial, 

FyFj7jFfe;F;i = 1. 

With these notions we can now formulate the key definition: 

Def. 2.4. Let f : M — > 5'^ be a discrete isothermic net with cross ratio factorizing function a. 
We say that 

(A, ij) ^ T^. := Ff7^°'^ e llom{{]} x R'^^\{i} x K^.i) (2.15) 



1,13) ^f f „,J,„^„ \ r- JO c,^^ i-Unt^ 1 -J- \,,.. 



defines the isothermic family of connections of f, where X € R so that 1 ^ Xaij for all edges 



Note that the F'*' are metric connections on M x R*'^, as all F^ are isometrics and hence descend 
to connections on M x S^. 



We alread y know that , if / is an isothermic net, then the isothermic family of connections 
is flat, see 



2.15 



18, §5.7.5]. Here we shall give an independent proof, not relying on the quaternionic 



setup, as well as a certain converse of this fact (cf. [17, Thm. 3.14] ) 



11) Wc shall make the notion of a fiat (discrete) connection precise in the following section. 

12) Here M6b(3) acts on S"^ = Im JET U {oo} by Mobius transformations, i.e., by fractional linear transformations. 

13) We have T^T^j = 1 on all edges (ij), so that the F^ qualify as discrete (linear) connections. 
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Lemma 2.5. Let f : M ^ S^ be a regular discrete net, i.e., any three vertices of a face are in 



general position, and let a be a function on the edges. Then the connection given by (2.15] is flat 
if and only if f is isothermic with cross ratio factorizing function a. 

Proof. First note that F,. ,°'^(X) = X mod /^ when X J. fi, so that F,. ,°'^ projects to the 
identity on fj^/ft and similarly for fj. Consequently, 

The same is true for the product Ff~,°''F j^,""" so that, if we now assume flatness of the connec- 
tion, 

fiJi fjJk fijl fljk ' ' 

we learn that F"^ — id on [fj ® fi)^ since fj ^ fi. Moreover, fj and /; are eigendirections of F'^ 
and 

■X ax e 



r^W = <! X iix±f„fi, 

^r^X iiXefi. 

l — Xau ■' ^ 

As F^ is, along with F , j,"'^ and F , ^ ^'^ , an orientation preserving orthogonal transformation we 
infer that, for all A, 

i-Aa,, i_A,,, ^ ^ ^ F^ = Fy''(^\ g(A):=ii^. 

l-Aoij l-Aoii JjJi ' ^V J 1-Aajfc 

Hence a^ = a^fe and a-a — aik (in which case F"^ = id for all A) or Uij = aik and ajk = an. 
Now decompose Fk £ fk\ {0} as Fk = Fi + Fj + F,}- e /, ® fj (/i ® J^)^ and observe that 

pA/'p "i 1 j^l — Aajj ^ ^ N 1 — Aajj p _^ 1 p_L j 1 p °-ij p 

i l^fej - l-\a^k^ fijj ^^''' ~ 1-Aa,fc ^« ^" l-Aa^fc -^fc "I" (1-Aa,j )(1-Aajfc ) -^J a^^ -^» 

as A -^^ cxD. This shows that, since fk j^ fi, we cannot have a^ = ajk] hence the function a has 
equal values on opposite edges of an elementary quadrilateral. Moreover, we learn that 

h = rf:}!if.\ 

showing that the four vertices fi, fj, fk and fi are concircular and the edge function a factorizes 
their cross ratio, 

[fi-Jj'Jk-Ji] =q{oo) = ^. 

J"- 

Hence / is a discrete isothermic net. 

Conversely, suppose that / is discrete isothermic with cross ratio factorizing function a; we wish 

to show that 

pA pA _ pl/?(^) _ j^qW _ pApA 

^ y^ jk - ^ fjji - ^ fijj - ^ a^ Ik' 
where q{X) ~ -^~^°'^ . As the second equation holds and the third is obtained from the first by 
exchanging the roles of j and / it suffices to prove the first of these equations. Also, 

p9(A) pA pA / x.^ _ / ^ if X e fj 

^ f, Ji ^ *i^ 3k{^)- <^x mod /j ax ± fj 

so that flatness of the family of isothermic connections of / follows as soon as F*^ I T^F^j, has 
another isotropic eigendirection — we shall show that fk serves this purpose: consider 

jRpi = JRU {oo} 9 A .. rf^r^kifk) = r,i^(/fe),r}/y^(/fe) e s\ 

Both maps parametrize the same circle in terms of a certain cross ratio, given as a linear fractional 
transformation of A; in particular, 

A = ^ T^^ifk)^Tl^^{fk)^fk^rljSfk)=T](':i^{fk), 
A = oo ^ F^^. (/,) ^ F- ,^. ifk) = /. = F,f,, (/,) = F}/f ) (/,), 

^ - ^ ^ r^(A) - r?.,,/AO = /, ^ r-j,{h) = T)/f\fk). 

As two Mobius transformations of a circle coincide as soon as they coincide at three points we 
conclude that F^^. (fk) = F}/ J^^^ (fk) for all A. < 

Note that, freeing the second part of the proof from the specific notations of the situation, we 
have proved the following: 
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Lemma 2.6. Write [pi; P2', P3', Pa] — f with a,b £ R for the cross ratio of four concircular points 
Pi e S^, i = 1,. . . , 4; then, for all A e JR, 

pl-aApl-6A _ p(l-bA)/(l-aA) _ -p(l-aA)/(l-6A) _ yl-bXyl-aX 

Pl;P2 P2;P3 -*- P2 ,P4 -•- P4 ,P2 '^ pi ,pi '' pi ,p3 ■ 

Thus, for a discrete isothemiie net /, there are gauge transformations T^ : M -^ 0(4, 1) identifying 
the connections F'^ on M x IR^'^ with the trivial connection: 

Lemma & Def. 2.7. Let f : M ^ S^ be a discrete isothermic net with its isothermic family of 
connections T^. Then the gauge transformations 

T^ : M -^ 0(4, 1) with T/ = T^^T^^ 

are the Calapso transformations of f; the isothermic nets f^ := T^ f are its Calapso transforms. 

Note that a^ = jz^ is a cross ratio factorizing function for the Calapso transform /^ of / with 
cross ratio factorizing function a, sec [18, §5.7.16]|; the isothermic family of connections of /^ is 



given by 

- ij — ^i ^ ij [-'- j 



TtTf+\Tp-\ (2.16) 

which shows that the Calapso transformations of a discrete isothermic net satisfy a 1-parameter 



group property, see [17] or [18, §5.7.30] 



rpfl,Xrpfl _ y/J + A (2.17) 



3. Polynomial conserved quantities 

The second key notion in our definition of discrete cmc nets in space forms will be that of poly- 
nomial conserved quantities: 

Def. 3.1. Let f : M ^ S^ be an isothermic net. A polynomial conserved quantity of f is a map 

RxMb (A, i) ^ P,(A) = ^;Lo /^("U^- e ffi4.i[A] 

so that, for every fixed A, 

T^P(X) = const. 

Hence, a polynomial conserved quantity of an isothermic net can be thought of as a polynomial 
family of parallel sections of the vector bundle M x IR^'^ equipped with the isothermic family of 
connections: 

T^P{X)= const. ^ P,{X)=r^jPj{X) (3.1) 

on all edges (ij) of M. 

3.1. Basic properties 

Clearly, as T"^ acts linearly on IR"*'^, the polynomial conserved quantities of a given discrete 
isothermic net / can be superposed: 

Lemma 3.2. The space of polynomial conserved quantities of f is a vector space. 

As a consequence, we can construct new polynomial conserved quantities from a given one by mul- 
tiplying with real polynomials p(A), thereby raising the degree; for example, if P is a polynomial 
conserved quantity of / then (1 + A)P(A) will be a new polynomial conserved quantity of higher 
degree. Thus we will be interested in (non-vanishing) polynomial conserved quantities of lowest 
possible degree. The following lemma provides a criterion: 

10 
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Lemma 3.3. Let P{^) = for a polynomial conserved quantity P{X) : M — ^ JR"''^[A] of f; then 

P(A):=^P(A) 
is a polynomial conserved quantity of f of lower degree. 
Note that, if Pi{fi) = for some i e M then P(^) = on M smce T^'P{^i) = const. 

Proof. Writing Pi{\) £ K^'^[A] in terms of a basis of IR'^^^ shows that /i is a common zero for all 
(real) component polynomials, which are therefore divisible by (A — /i) . Hence Pi (A) is polynomial 
at any i e M . 

Clearly 

r^P(A) = ^ T^P{\) = const 

for any fixed A, showing that P(A) is a polynomial conserved quantity of /. < 



As a direct consequence of the previous two lemmas we learn that, if two distinct polynomial 
conserved quantities P{\) and P(A) of degree A^ G W of an isothermic net have the same value at 
some point (/i, i) £ Rx M, then there is a polynomial conserved quantity of degree < TV — 1. In 
particular: 

Cor. 3.4. A non-zero polynomial conserved quantity of lowest possible degree N , 

P(A) = X^Z + . . . + A"g : M ^ K^'MA], 
is uniquely determined by either its top or bottom coefficient Z or Q, respectively. 

Since the T^ are orthogonal transformations, there is another obvious property of a polynomial 
conserved quantity which will become important later: 

Lemma 3.5. If P(A) is a polynomial conserved quantity of f, then |P(A)p depends only on A; 
in particular, \Z\'^ and jQp are constants. If P{X) = XZ + Q is a linear conserved quantity, then 
also {Z, Q) = const. 



Finally note that the equation (3.1) for a polynomial conserved quantity depends crucially on the 
choice of a cross ratio factorizing function a: however, if 5 := aa is a new cross ratio factorizing 
function, then 



by (2.15); hence P(A) = P(aA) is a new polynomial conserved quantity satisfying (3.1) with the 



new isothermic family of connections. Consequently: 

Lemma 3.6. If P{X) is a polynomial conserved quantity of f with respect to a as a cross ratio 
factorizing function then 

P(A) := P(aA) 

is a polynomial conserved quantity of f with respect to a := aa as a new cross ratio factorizing 
function. 

3.2. Geometric properties 



We now turn to a more detailed analysis of (3.1)| and its geometric consequences: using (2.4) the 



fact that a polynomial conserved quantity P(A) is a family of parallel sections of the isothermic 



family of connections (2.15) of / reads 

P,:(A) = r^^.P,(A) = P^.(A) + ^{^4^(P,(A),P.)P,-(P,(A),P,-)PJ 

on any edge (ij) of M; exchanging the roles of the endpoints i and j of the edge we obtain a 
similar equation which, as fi ^ fj, yields two equations 



dP^,W = j^{{P,{X),F,)F,-{P,{X),F,)F,} 



1- 



for some light cone lift F of /. Note that the second equality follows from the first by taking 
scalar products with Fi and Fj, respectively. Hence, we also obtain the converse: 

11 
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Lemma 3.7. P{X) is a polynomial conserved quantity of f if and only if, for all edges (ij) in M , 



Note that, in case F is a Moutard lift of / satisfying 1(2.5) , then (3.2)| simphfies to 

dP,j{X) ^X{p,{X)F,^p,{X)Fj}, where p{X) := (P(A),F>. 
The integrabihty d^P{X) = of this equation then yields 

{pk{X) - p^{X)){F, - Fi) = {p,{X) - pi{X)){Fk - FO, 



(3.2) 



and hence p(A) satisfies the very same Moutard equation (2.6) as F does 



Now the key observation from (3.2) is that this equates a polynomial of degree N and a polynomial 
of degree A^ + 1 with vanishing constant coefficient. Hence, looking at the degree and degree N 
and A^ + 1 terms we obtain the following two corollaries: 



Cor. 3.8. IfP{X) ^X^Z 



Q is a polynomial conserved quantity of f, then Q — const. 



Thus a polynomial conserved quantity natu rally provides a n ambient quadric Q of constant cur- 
vature K = — IQp for the isothermic net, see [18, Sect. 1.4] : 



Q^{YeL^\{Y,Q)^-l}. 



(3.3) 



Cor. 3.9. If P{X) = X^ Z + A^ ^F + . . . + Q is a polynomial conserved quantity of f then: 

(i) Zi _L Ji at all points i in M; 

(ii) Zi + aij (f^j') Fi = Zj + ay -j^^ Fj for all edges (ij) in M; 
(iii) |Z|2 > and \Z\^ ^ if and only if Z e /. 

Here, (iii) follows from (i) since /^ carries a positive semi-dcfinitc metric with only / as a null 
direction. Also note that, since |Zp is constant, cither Z || F at all points or, without loss of 
generality, \Z\'^ = 1 as the space of polynomial conserved quantities is hnear. 

We will be mostly interested in the latter case. To interpret this situation geometrically first note 
that, if |Zp = 1, then i i— s- Z^ defines a discrete sphere congruence so that every sphere Zi contains 
the point fi, by (i); moreover, the sphere 



^ij '• — '^i 



,14) 



"*J {Fi,Fj) ^'^ ^ ^J ^ "U {Fi,Fj) -^ 



(3.4) 



belongs to both contact elements ' defined by Z at the endpoints of an edge so that 

M5t^Z, + f^:= {Z, + aF,\aeR} 



defines a discrete principal net in Lie geometry, see [6, Sect. 4.1] , with curvature spheres Sij . On 
the other hand, the existence of the curvature spheres Sij, which touch both spheres Zi and Zj, 
can be interpreted as a discrete version of the enveloping condition for a sphere congruence defined 
at the vertices of a discrete net. Thinking of / as a net in JR'^, the spheres Z define a unit normal 



field n at the vertices of /, which satisfies the trapezoid property of [22, Sect. 3] 
We summarize these observations in the following ' 



14) Recall that, in contrast to Mobius geometry, Lie geometry considers oriented spheres; thus ±Z + / will define two 
contact elements which differ by the orientation of their spheres. 

15) We insist on a consistent orientation of the spheres of an enveloped sphere congruence. 



12 
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Cor. &: Def. 3.10. If P{X) = X^ Z + . . . + Q is a polynomial conserved quantity of f with 
|Zp = 1, then f envelops the discrete sphere congruence Z: we say that f : M ~* S^ envelops a 
discrete sphere congruence S : M ^ S^'^ if 

(i) Si e fj^ for all i e M (incidence) and 

(ii) Sj = Si mod fi ® fj for each edge (ij) of M (touching); 

the common sphere Sij of the two contact elements Si + fi and Sj + fj given by Si and Sj at the 
endpoints of an edge will be called a curvature sphere. 

Note that the condition for a sphere congruence to be enveloped by a net is a condition on the 
congruence of contact elements defined by the sphere congruence — hence, if / envelops S and F 
is any light cone lift of /, then any sphere congruence S + hF, h : M ^ M, is also enveloped by /. 

Before focusing on this geometrically interesting configuration we shall, for the rest of this section, 
investigate the degenerate case where we can, without loss of generality, take Z — F as a canonical 
light cone lift of /. 

Lemma 3.11. If P{X) = X^ F + . . . + Q is a polynomial conserved quantity of f, then F is a 
Moutardlift off. 



Proof. From (ii) of Cor. 3.9| we learn that 



{F„Fj) + a,j{Yj,F,) = {F„F,) + a,,{Y„F,) = 



since Fi and Fj are linearly independent; thus ' {Y, F) = const ==: c and F satisfies (2.5)| with 
the cross ratio factorizing function a :— ~ca. < 

As an example we seek an isothermic net with a degenerate degree 1 polynomial conserved quan- 
tity: this is the lowest degree possible since the top term of a degree polynomial conserved 
quantity is constant and can therefore not be a (Moutard) lift of an isothermic net. 

Now consider, as before, 

|(xo,...,a;4)|2 = -a;2^^ti^f 

as the quadratic form of the Minkowski scalar product of JR^'^ and the isothermic net 

{-1, 0, 1}2 3 (m, n) ^ /(„,„) := (77m, i±^ + (-l)"ifi, /?n) G R^ (3.5) 

where a S (0, 1) and /?, 77 > 0, and let 

^(™,„) :- (-l)"i^(™,„) and Q := -^(l±s, 0, 1, 0, -i±2), (3.6) 

where F = (^-^5^, /, ^— 5^) is a Euclidean \iit^^\ Note that IQp > so that Q describes a 

sphere and / appears to be a perturbation of a net on that sphere. 

Since 



and / has rectangular faces, so that the cross ratio (2.3) on an elementary quadrilateral becomes 



[jii fj: fk', fl\ — l/i-/, P ' 



16) We neglect the case where Z may have zeroes, i.e., where the degree of P(A) may not be constant. 

17) This we already knew from Lemma 3.5, because {Y,F) is the A^^~^-coefEcient of |P(A)p. 



18) Cf. 1(2. 7j and |(2.9)| : this is the Euclidean lift with respect to Qg = (1, 0, 0, 0,-1), defining a fiat quadric of constant 
curvature via |(3.3j[ 

19) Suppose that |Qp < for a degenerate linear conserved quantity P{\) = \F + Q of an isothermic net /; then 

aij = {F„Fj) = -^\dFij\'^ <0 



for any edge {ij) since dFij ± Q by Lemma 3.5 and / is regular. Hence the cross ratio of any face becomes positive; thus, 



if we seek an isothermic net with embedded faces, then Q necessarily describes a sphere. 

13 
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Oy := \{Zi, Zj) provides a cross ratio factorizing function, i.e., Z is a Moutard lift of /. Moreover 

{XZ + Q,Z) = {Q,Z) = -2, 
so that 



d{XZ + g)y - j^^{{XZ + Q, Z)jZ, - {\Z + Q, Z),Zj) = 



and P(A) := \Z + Q is a linear conserved quantity for / by Lemma 3.7, 
Note that Z = P(oo) is a lift of / and that 

is a lift of a Mobius equivalent net or, more precisely, of the "antipodal" net in the quadric of 



constant curvature given by Q, see [18, Section 1.4] . We shall come back to this observation later. 
Fi nally o bserve that our restriction of the dom ain t o {—1,0, l}'^ was not necessary: /, as given 



by (3.5) , can be extended to all of Z^ while (3.6)| keeps defining a linear conserved quantity 



for /. However, this restriction will be convenient later when we shall recycle this example to 
demonstrate another aspect of our theory. 

3.3. Special isothermic nets 



As in the smooth case, see [21, Sect. 2.2] , we can now use the existence of a polynomial conserved 
quantity to define a special class of discrete isothermic nets, ordered by the (minimal) degree of 
an associated polynomial conserved quantity. However, in contrast to the smooth case, where 
a polynomial conserved quantity is essentially unique because its top degree coefficient encodes 
the conformal Gauss map of the underlying isothermic surface, the space of polynomial conserved 
quantities of minimal degree may be higher dimensional ' and may contain elements with null top 
degree coefficient, as in the above example, even though we require the existence of a normalized 
polynomial conserved quantity of minimal degree: 

Def. 3.12. A polynomal conserved quantity P{\) — \^ Z + . . . + Q of an isothermic net f will 
be called normalized ii \Z\'^ = 1; we say that f is a special isothermic net of type N if it has a 
normalized polynomial conserved quantity of degree N, but not of any lower degree. 



As a direct consequence of this definition and the 1-parameter group property |(2. 17) of the Calapso 



transformations we obtain stability of the class of special isothermic nets of a fixed type N under 
the Calapso transformation: 

Thm. 3.13. If f is special isothermic of type N then so are its Calapso transforms f^ — T^f. 



Proof. T/^'^ = yM+A(7-M)-i are the Calapso transformations of f by |(2.17)| , see also [18, §5.7.30] 



Hence, if P(A) = X^ Z + . . . + Q is a polynomial conserved quantity of /, then 

P^(A) :=r^P(Ai + A) 
defines a polynomial conserved quantity of /^ of the same degree; moreover, 

|P'^(A)|2 = |P(A + m)|2 = A2^|Z|2 + A2~-i ... + ... 

showing that P^(A) is normalized as soon as P(A) is. Thus f^ is special isothermic of type < TV. 

On the other hand / — T^~^^f^^ is a Calapso transform of /^, by |(2.17) again, so that the same 
argument shows that /^ is special isothermic of type > A^. < 

So far it is rather unclear how restrictive the definition of special isothermic of a given type is: as 
we shall see later (and as one might expect after an equation count), the condition on a discrete 
isothermic net of being special isothermic of type A^ is not a local condition — in particular, we 
will see that every isothermic 3 x 3-net has plenty of linear conserved quantities and is, therefore, 
special isothermic of type 1. On the other hand, the condition of being special isothermic of type 
does already impose a condition on a 3 x 3-grid, and the corresponding constant conserved 
quantity is generically unique: 



20) However, we do expect uniqueness of a (normalized) polynomial conserved quantity of minimal degree for generic 
isothermic nets of a sufficient size. 
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Thm. 3.14. / is special isothermic of type if and only if f takes values in a 2-sphere. 

Proof. First suppose that / is special isothermic of type with (necessarily ' normalized) constant 
conserved quantity P{X) ~ Q. T hen Q defines a fixed 2-sphere, since |Qp = 1, that contains the 



points of the net by Cor. 3.9 (i) . This also shows that a (normalized) constant conserved quantity 
is unique as soon as the isothermic net does not take values in a circle. 

Conversely, if / takes values in a fixed 2-sphere S C S^ then the Calapso transformations T^ 
of / take (up to a constant Mobius transformation) values in the Mobius group Mdb{S) of this 



2-sphere (see [17, Thm. 3.14] or [18, §5.7.22] ). That is, the Calapso transformations T of / can be 
chosen to fix a unit vector Q defining the 2-sphere S: this yields a normalized constant conserved 
quantity so that / is special of type 0. < 

As a more involved example we take up surfaces of revolution and investigate the symmetry 
of a corresponding polynomial conserved quantity: suppose F(m,„) = (— l)'"(Afm + $„C) is the 



Moutard lift of a discrete surface of revolution, see (2.8), with a rotationally symmetric polynomial 
conserved quantity, i.e., we assume that $~^P(m„)(A) ~: P„i{X) does not depend on n. Then, 
clearly, 

P(m,„)(A) := (P(™,„)(A),F(„,„)) = (-l)™(P„(A),A/™ + C) 

is independent of n; we shall see that the converse also holds: 

Lemma 3.15. A polynomial conserved quantity P{X) of a discrete net f of revolution with 
canonical lift -F'(m,n) = (— l)™^ri(Afm + C) is rotationally symmetric, P{rn.n){^) = ^nPm,{X), if 
and only if 

p(A):-(P(A),F) 

does not depend on n. 

Proof. We already know that p{X) depends only on m if P{X) is rotationally symmetric; to prove 
the converse we assume that p{X) does not depend on n and write 

P(m,n) = -P(i,„) + <^n{a{m.n)C + /3(„,„)C^) E R^'^ R^ , (3.7) 

where C'^ complements C to form an orthonormal basis of R^ and a and (3 are suitable real valued 
functions. First note that 

(~1)>„ = (^(i,„), Af,„) + a(„^„). 

Now we fix m and, in order to simplify notation, consider all functions as functions of n only. 
From (3.2) we get 

= dPn,n+l{X) + Xp{X)dFri^„+i; 

hence, in particular, 

d'Pn,n+i = a-iid dan,n+i = -(d/',^„+i,M) = 0. 

For a function g on the vertices let gn,n+i '■= ^" 2"^^ denote the associated function on the edges, 
and let dipn,n+i be the rotation angle along the edge, i.e., the rotation angle of <i>„+i$^^; observe 
that 

d(*C^)„,„+i = -2tan%^($C)„,„+i ± dP„,„+i, 



2($C^)„,„+l = cot ^i^^^ d($C)„,„+i II dFn,n+l- 

Hence the JR^-part of = dP„^„+i(A) + Xp{X) dFn^n+i yields 

= /3„,„+i(A) and = a(A) + (-l)"Ap(A) + id/3„,„+i(A) cot ^ 



Thus, considering two consecutive edges, = /?„ sin "^'^ ^'"^ '^"■"+'- so that j3 vanishes and 

P(„,„)(A) = P,^(A) + (-l)"+iAp„(A) $„C (3.8) 

is clearly rotationally symmetric. < 

As a simple consequence we see that a linear conserved quantity of a discrete surface of revolution 
/ is rotationally symmetric as soon as / is rotationally symmetric in the space form defined by 
the constant term: 



21) Remember that a constant conserved quantity cannot be degenerate by Cor. 3.9 (i), since / is not constant. 
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Cor. 3.16. A linear conserved quantity P{X) = XZ + Q of a net f of revolution is rotationally 
symmetric if and only if Q is. 



Proof. By (3.8), the constant term of a rotationally symmetric polynomial conserved quantity has 
no JR^-component in the decomposition (3.7); for the converse observe that, in the case of a linear 
conserved quantity, P(,„,„) = (Q, i^(m.n))- < 



4. The Backlund transformation 

Bianchi's Backlund transformation of smooth constant mean curvature surfaces in Euclidean space 
turned out to be a special case of the Darboux transformation of isothermic surfaces: considering 
the Darboux transformation as an initial value problem depending on a real (spectral) parameter, 
the Darboux transforms of a constant mean curvature surface turn out to have constant mean 
curvature as soon as a certain relation between the initial value and the parameter is satisfied, see 
18, §5.4.151; in fac t, they turn out to be the Backlund transforms of the surface, see [15, Thm. 7] 



sense of |[16] 



and [20, Thm. 4.4] , A similar fact holds true for discrete constant mean curvature surfaces in the 



In [21, Thm. 3.2], this condition for a Darboux transform of a constant mean curvature surface to 
have (the same) constant mean curvature again is shown to be a special case of a similar condition 
for (smooth) special isothermic surfaces of given type N. Here we shall analyze the situation for 
discrete special isothermic nets, giving rise to what we will call the "Backlund transformation" 
for special isothermic nets. 



To this end we shall first recall the Darboux transformation for discrete isothermic nets, cf. [16 
or 



[18, §§5.7.12 fc 5.7.19] 



Def. 4.1. Let f : M ^ S^ be a discrete isothermic net with its family of Calapso transformations 
T^ . Then a Darboux transform of f is a discrete net 

f : M -^ S^ so that 3fi e R : Tf^f = const. 
Equivalently, we can characterize a Darboux transform by the condition 

f "pM f pi — aij/^ r 



for all edges (ij) of M: remember that the isothermic family of connections (2.15) descends to 
a family of connections on M x S^; thus a Darboux transform / can be thought of as a parallel 
section of a connection F'^ on M x S^ of the isothermic family of connections. Using a cross 
ratio identity (see [18, §4.9.11 ), this condition can be reformulated as the cross ratio condition ' 



(cf. (2.4)) 



[Jii Jji Jj'j Ji\ ^ 1 [/ji /ij Ji'^ Jj\ ~ 0,ij^. 



(4.1) 



Clearly, a Darboux transform / of an isothermic net has a light cone lift F : M ~* L* so that 

T'^F EE const. ^ F,^ T%Fj (4.2) 



on all edges {ij) of M, t hat is, F is a F^-parallel lift of /. Note that this is exactly the conserved 
quantity condition (3.1) for a light cone map F and a fixed value /i of the spectral parameter; we 
shall come back to this point later. 

The crucial observation is that the Darboux transformation produces discrete isothermic nets 
from isother mic nets, see [16] or [18, §5.7.12], where a proof relying on the he xahedron le mma 



[18, §4.9.13] was given; here we shall again give an alternative proof, relying on Lemma 2.5 , that 
als o pr ovi des us with a useful formula for the Calapso transformations of a Darboux transform, 
cf. [it] or |[18, §5.7.35]| : 



22) Using quaternions or the Clifford algebra of JR^ and thinking of 5^ = IB? U {cxd}, this is equivalent to the (discrete) 



Riccati type equation 



dh=i^U-f),df;Af-fh 



where df*- = aij(dfij) ^ is the derivative of the Christoffel transform of /, cf. [18, §5.7.7] ; the condition fi = T''^A'j appears 



as the usual linearization from this Riccati equation: as "Darboux's linear system" . 
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Lemma 4.2. A Darboux transform f, T'^f = const, of a discrete isothermic net f with cross 
ratio factorizing function a is isothermic with the same cross ratio factorizing function a = a and 
with Calapso transformations 

(4.3) 



'T^A T^Apl -^/M 



/,/ 



Proof. We shall show that the isothermic family of connections |(2.15) of / is given by 



r^ := rn-^ = (r;^/-)rir^(r;;^/-)„ 






which is then, with F^, clearly flat so that / is isothermic with cross ratio factorizing function 
a = a by Lemma 2T^ ; further, the formula for the Calapso transformations T'^ also follows directly. 

Thus we have to show that 

-pi — A//^-pl — Aaij -pi — Aa^j -pi — A//^ 



' f^J^ f,Jj 



' ii'fi fjjj 



But, with (4.1) , this follows directly from Lemma 2.6 . < 

At this point we are now in perfect shape to discuss: 

4.1. Darboux transforms of special isothermic nets 

As outlined at the beginning of this section, we aim to obtain the Backlund transformations for 
special isothermic nets of type N as particular classes of Darboux transformations. The following 
theorem provides the essential criterion: 

Thm. 4.3. Let f be a Darboux transform, T^f = const, of a special isothermic net f of type N 
with normahzed polynomial conserved quantity P(A) of degree N . Then: 

(i) / is special isothermic of type < A^ + 1 and 

(ii) / is special isothermic of type < N as soon as P{fJ.) G /^. 

Note that, if Pi(/i) G f.j^ for some i £ M, then this holds true for all i G M. Namely, if i^ is a 



light cone lift of / as in |(4.2)| , then P{ii) and F are both parallel sections of the metric connection 



r^- on M X IV'^ so that their scalar product {P{fi),F) = const. 



Proof. We define 



P(A):=(A-/i)r)r//'^P(A). 



(4.4) 



Clearly T P(A) = const by (4.3) and, using (2.4), we obtain 



P(A) = (A - A^)P(A) - ^{ APi_ii) (p(A), F)P + A(P(A), F)F}, 



which is polynomial of degree < A^ + 1 as (P(A), F) has degree < A^ — 1 by Cor. 3.9 (i). Finally, 
writing 

P{X) = X^Z+... + Q and P(A) = A^+^Z + . . . + Q, 



we find that 



A2^+2|zr 



= |P(A)|2 = (A - ^I)^ |P(A)|2 = A2^+2|Z|2 + 



for all A, so that |Zp = \Z\'^. Hence P(A) is a normalized polynomial conserved quantity of degree 
A^ + 1 and, therefore, / is special isothermic of type < A^ + 1, proving (i). 
To prove (ii) note that 



when P{fJ,) G /^. Hence / has a polynomial conserved quantity of degree < Af by Lemma 3.3; in 
particular, 

P(A) := Fjr^/'^P(A) (4.5) 
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provides a normalized polynomial conserved quantity of degree N for / showing that / is special 
isothermic of type < A^. < 

Def. 4.4. A Backlund transform f of a special isothermic net f of type N with polynomial 
conserved quantity P{X) is a Darboux transform, that is, T^^f = const, so that P{ijl) G /^. 

The Backlund transformation between special isothermic nets of the same type is symmetric, 
that is, if / is a type N Backlund transform of a special isothermic net / of type N with polynomial 
conserved quantity P{\) then / is a Backlund transform of /, where P{X) is given by (4.5) : firstly, 
/ is a Darboux transform of / since the cross ratio condition (4.1) is symmetric in / and /, 



VJ^'^ Jj'j Jj\ Ji\ — Uii Jjj Jj\ Ji\ — 0,ij^, 



by a cross ratio identity, see [18, §4.9.11] ; and secondly, using (2.4)| , we find that / satisfies the 

F mod /. 



Backlund condition^'*' -P(m) G /^: 

P{fi) = hm^^^ P(A) = P{fi) - ^^i^^ 



(^,^> 



Recall that constructing a Darboux transform / of a given isothermic net amounts to determining 
a parallel isotro pic section F of M x R*'^ equipped with a connection T^ in the isothermic 
family, see (4.2) : after choosing an initial value there is a unique solution. As the condition on a 
Darboux transform of a special isothermic net to become a Backlund transform is preserved by 
the propagation it is sufficient to choose the initial value of the Darboux transform appropriately 
in order to obtain a Backlund transform. 

Now note that (4.5)| has exactly the same structure ' as the con served quan tity condition (3.1) , 
with - taking the role of the a^. Thus, completely analogous to Lemma 3.7 , 



pw + ^,^W^^ = pw + ^,W^p^ 



(4.6) 



26) 



where F and F are any ' light cone lifts of / and /, respectively. In particular, writing 



P(A) =X^Z + X 



\N-1 



Y 



Q and P(A) = A^Z + A^-iy 



Q, 



we find that, analogous to |Cor. 3.- 



geometric interpretation of Cor. 3.9 



the constant coefficients Q = Q and that, analogous to the 



S 



JY.F) 
t^iF.F) 



F = Z 



(Y.F) 
l^(F.F) 



F 



(4.7) 



27"! " 1^ I 

yields a sphere congruence ' that is enveloped by both / and / in the sense of pef. 3.10| . Thus our 
Backlund transformation for special isothermic nets is a sp ecial type of Ribaucour transformation 
for discrete curvature line nets in Lie sphere geometry, see [6, Def. 22] . 



23) Wc know that the Backlund docs not increase the type but, in special circunistances, if may decrease the type (by 1). 
We shall come back to this point later. 

24) The assumption that / is special isothermic of the same type N as f ensures that P{\) has, with P(A), minimal 
degree and, in particular, that -P(^) ^ so that the condition is meaningful. 



25 
cf.|[6] 



f'his is a common phenomenon in discrete differential geometry, often referred to as "multidimensional consistency" , 
or |7|: the transformations of a discrete net are governed by exactly the same conditions as the net itself 

26) A canonical choice would be to take F and F to satisfy (2.5) and {F, F) = j^'- this is the "3D-consistency" of the 
condition to be discrete isothermic, see [7l 

27) Indeed, a similar fact can be proved for Darboux transforms in general: given a sphere congruence Z enveloped by 
an isothermic net and a Darboux transform / of /, the sphere congruence 



S:=Z- 



{Z,F} 
{F,F} 



will be enveloped by both f and /. This is a fact about the existence of Ribaucour transforms for discrete curvature line 
nets in Lie geometry, cf. r6]l Note that this is a different "Darboirx sphere congruence" than the one discussed in [16] , 
which "lives" on the faces of the domain M. 
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Thm. 4.5. Let f be a Backlund transform of a special isothermic net f of type N, with respective 
polynomial conserved quantities 



PiX) ^ X^ Z + X^-^Y 

Then 
(i) Q = Q, and 

(ii) Z + 



_ pl-A/^ 



and P{X) = r'-^'^PiX) = X^Z + X^-^Y + ... + Q 



{Y,F) p^2 ' <^-P^> 



F; and, 



f,{F,F) f,{F,F) 

iii) in particular, Z + f and Z + f give rise to a Ribaucour pair in the Lie geometric sense. 



4.2. Bianchi permutability 

A key feature of Darboux-Backlund type transformations for smooth or discrete classes of sur- 
faces is Bianchi permutabihty: given two transforms of a surface there is a fourth (often unique) 
surface, which is a simultaneous transform of the two initial transforms — thus providing the 
combinatorics of a quadrilateral for the transformation. We will refer to such quadrilaterals as 
Bianchi quadrilaterals. 

In particular, such a theorem holds true for the Darboux transformation of (discrete) isothermic 
surfaces, see [16] or [18, §5.7.28], where the fourth surface is uniquely determined ' by a cross 
ratio condition : 



Thm. 4.6. If fi and /2 are two Darboux transforms, T^^ fi = const, of a discrete isothermic net 
f, then 

~.fJ.2/l-ll 
' fuf2 



h 



pM2/A'l £ 



is a simultaneous Darboux transform of /i and /2 .' 

Ti^ fi2 = const and r^^/12 = const. 



We wish to prove a similar theorem for the Backlund transformation of special isothermic nets of 
a given type N. 

Thus let / be a special isothermic net with polynomial conserved quantity P(A) and let fi be two 
Backlund transforms of /, i.e., 

T^'/,; = const and P(^,) € /^. 



Clearly, by the permutability theorem Thm. 4.6 for the Darboux transformation and (4.5) , the 
simultaneous Backlund transform /12 of /i and /2 is necessarily 



_ -pfJ.2/fJ.l 



_ pl-A/^2 



_ -pl-A//ii 



/i2 = r^^'Tf with Pi2(A) = r;-^ ''^Pi(A) = r'-^^^p^ix) 

J1,J2 J12,J1 J12,J2 



(4.8) 



as a polynomial conserved quantity — just as in the proof of [Lemma 4.2| it follows from [Lemma 2.6 
that 



pl-A//i2pl-A//ii 

/l2,/l /l,/ 



r 



(1-A/mi)/(1-A/m2) 

' fl.f2 



pl-A//ji-pl-A//j2 

fl2j2 /2,/ 



SO that Pi2(A) is well defined by (4.8) . Now it is straightforward to see that 



Pi2(mi) = limA^A^i ^i f. 



(1-A//.i)/(1-A/m2) 



p(A)^p(m)-^^;^^A 



■he fir, 



showing that /2 is a Backlund transform of /12; similarly, /i is also a Backlund transform of /12. 
The symmetry of the Backlund transformation then completes the proof of the following 



28) This is in cont ast to the Ribaucour transformation of (discrete) principal nets, where a l-paramctcr family of fourth 
surfaces exists, see 



[10 or 5. 



29) Note how repeated application of this theorem builds up a discrete isothermic net. 



19 



discrete erac surfaces 17 April 2008 - ^9 arXiu uersiou 

Thm. 4.7 (Bianchi permutability). Given two Biicklund transforms /i and ji with parameters 
Hi and ii2, respectively, of a special isothermic net f of type N, there is a net /12 so that the four 
nets form a Bianchi quadrilateral: that is, /12 is a Biicklund transform of /i with parameter fj,2 
and of /2 with parameter /ii. 

Note that we have not used any new arguments to prove this theorem — indeed, using the simi- 
larity of the polynomial conserved quantity equations and the conditions governing the Backlund 
transformation, we could have formulated a proof based on the fact that "3D-consistency" of a 



■'2D-systeni" implies higher dimensional consistency, see [7, Thm. 7] : in the case at hand we were 



interested in "4D-consistency" . Thus, any higher dimensional permutability theorems can now be 
proved by purely combinatorial arguments . For example, we can now argue that a "Bianchi 
cube" can be (uniquely) constructed from a special isothermic net and three Backlund transforms: 
the existence of the eighth Backlund transform is ensured by the very same fact that ensured the 
existence of a Darboux transform and the compatibility of the Backlund transformation with the 
construction, as discussed above. 

4.3. Complementary nets 



As we already noticed earlier, the equation (4.2) on a light cone map F to provide a Darboux 



transform of an isothermic net / is exactly the conserved quantity equation (3.1) for a fixed 
parameter /i. Consequently, any zero fi of |P(A)p provides a (light cone lift of a) Darboux 
transform F — P{fJ,) of an isothermic net / with polynomial conserved quantity P(A). Moreover, 
since (F,P(/i)) = |P(/x)p = 0, this Darboux transform will, in fact, be a Backlund transform 
of /. In this section we shall discuss the role of these special Backlund transforms of a special 
isothermic net. 

Def. 4.8. Let f be a special isothermic net with polynomial conserved quantity P(A); those 
Backlund transforms f of f given by F — P{fJ,), where |P(/i)p = 0, are the complementary nets 
off. 

Clearly, a special isothermic net of type N has at most 2N complementary nets and, as |P(A)p 
is an even degree polynomial, there may be no (real) complementary nets: for example, let / be 
a type 1 special isothermic net with linear conserved quantity 

P{X)=\Z + Q and H:=-{Z,Q), k:^-\Q\^; (4.9) 

then the complementary nets of / are given by 



P{H ± VH^ + k) = {H± ^/m + k)Z + Q (4.10) 

so that the number of (real) complementary nets depends on the sign of H^ + k. 

On the other hand, if enough complementary nets are known, then the corresponding polynomial 
conserved quantity can be reconstructed: first we observe that, given A^ + 1 distinct parameter 
values A = /j,o , . . . , ^jln , 

PW = Eto ^(Mn) n™^„ ^-^. (4.11) 

where P{fJ.n) are F'^" -parallel and, with «„ :— Y\„^^n -u ' ^^^ leading coefficient 

Z^En=0'^nP{pn)ef^ 

of P(A) is a (constant) linear combination of the P(/i„). These are the assumptions that will allow 
us to reconstruct a polynomial conserved quantity from iV -I- 1 suitable Darboux transforms: 



30) This is in contrast with the Ribaucour transformation, where 3-diniensional perinutabihty, i.e., a "Bianclii cube" 
theorem, is the critical case as there is no uniqueness in the Bianchi quadrilateral. 



31) Note that P{H ± \/H^ + k) ^ since / is special isothermic of type 1 so that P{\) cannot have zeroes. 
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Lemma 4.9. Let F"^ , n ~ Q, . . . ,N, be T'^''''^ -parallel sections of M x JR'*'^ for pairwise distinct /i„ 
and suppose that 



for some constants a„ e K and all i E M. Then 



^N 



^(A):=E:=o«n^"n™^„(A-M™) 



(4.12) 



is a degree N pol ynomia l conserved quantity for f with top degree coefRcient Z; if\Zi\^ — 1 at 
some i € M , then (4.12) defines a normalized polynomial conserved quantity. 



Note that, if |Zp > 0, then / envelops the sphere congruence Z in the sense of Def. 3.1C, as 
desired: incidence is given by assumption and touching is guaranteed because the sections F" are 
pMn parallel so that 

dZ^J = E^=o «" dPp^ = mod /, © /, . 

This also implies directly that |Zp = const on Af . 

Proof. It remains to show that P(A) is indeed a polynomial conserved quantity, that is, we wish 
to show that 



T^P,{X) - P,{\) - dP,,iX) + ^0^{^-L-^(P^iX),F,)F, - {P,iX),F,)F,} 







for all A and each edge (ij) of M. Clearly, this equality holds true for A = /i„, n = 0, . . . , N, since 

SO that it holds for all A as soon as we know that r^-Pj(A) — Pi (A) is a degree TV polynomial. By 
definition Pi (A) and Pj{X) are degree N polynomials, and (Pj(A), Fj) is a degree A^ — 1 polynomial 
by incidence, Zj ± Fj. Moreover, 



so that 



l-ai,A\^ J 



= ^ nm=o(i - «y/^™) e!Lo «"(^r: ^«) 

1-3 

= 
Pj{X), Fi) is also a degree iV — 1 polynomial. Hence the claim follows. 



Note that, in Lemma 4.9, we did not require the P" to be isotropic sections of M x ]R ■ , that is, 
we did not require them to be lifts of Darboux transforms of /. 

Now suppose that /i is a simple zero of |P(A) P and let P = P{i-i) denote (a lift of) the corresponding 



complementary net of /. Then, from |(4.5) 



A (PW,F) 



P(\\ — P(\) A [^{^),-f ! p( \ A [^{Aj,^y 



A {P(\),P(p)) 
)> 



P 



such that 



since 






M 1 = ^.^ (P(A),P(m)) p 



2(P(p),F) m"A^P A-/J ^ 



X—fj, 



2iEmJim^(x-,) 



PW-PM 

X—fi 



0.|P'(m)P = 



as A ^ /i. Consequently, if /i is a simple zero of |P(A)p, then / is a complementary net of /, that 
is, the notion of complementary nets is symmetric. 

If, on the other hand, /z is a higher order zero of | P( A) | , then P{n) = so that / is special isothermic 
of lower type than / by Lemma 3.3, Indeed, all type lowering Backlund transformations arise in 
this way: suppose that / is a type A^ — 1 Backlund transform of a special isothermic net / of type 
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N — or, otherwise said, / is a Darboux transform of / which is not a Backlund transform. Then, 
their polynomial conserved quantities are related by [(4.4) : 



so that 



P(A) ^ (A - M)r;/'^^(A) = (A - M)(P(A) - ^%P F) - A<W> p 



spans /, which is therefore a complementary net of /. 

Note that /z is a higher order zero of |P(A)p = (A - iif\P{\)\'^. 

We summarize these results: 

Lemma 4.10. Let f be special isothermic of type N with polynomial conserved quantity -P(A). 
(i) If iJ, is a higher order zero of |P(A)p, then F = Pif^) defines a type N —1 Backlund transform 

off. 
(ii) If f is a type N — 1 Backlund transform of f, then f is a complementary net of f, f 3 P(/i) 

for some /x, where jjl is a higher order zero of |P(A)p. 

As a consequence of this lemma, a special isothermic net of type N can have at most N Backlund 
transforms of type A^ — 1 — generically, a Backlund transformation is between special isothermic 
nets of the same type, and is therefo re sym metric as discussed above. For example, consider a 



special isothermic n et / of type 1 : by K4.10)| , / is a Darboux transform of a type net, that is, of 
a spherical net (see Lemma 3.14 ), if and only if H^ + k, — Q. 

In the remainder of this section we shall discuss geometric properties of complementary nets of 
special isothermic nets of type 1 and 2. 

First consider a type 1 special isothermic net with two complementary nets, i.e., H^ + k > in 



(4.10) ; also, we assume k ^ 0, excluding the degenerate case, where one of the complementary 
nets becomes constant. Then 

i^±:=Z+^Q, where ^± ■.^H±y/m + k, (4.13) 

provides F'^* -parallel light cone lifts of the two complementary nets so that 

Thus the two nets are Mobius equivalent and, more precisely, they are "antipodal" in the quadric 



of constant curvature given by Q, see [18, Section 1.4]. Moreover, the "orthogonal circles" 



6=^ :=span{F,Z,F±} 



of the Ribaucour pair ' (/, /*), see Thm. 4.5| coincide in corresponding points: 



cf = spanji^j, Zi, Fj^} = spanji^^, Z,, Q}. 

For a type 2 special isothermic net / we obtain two similar properties: first let 

F'' ^ Pifir.) = fJ^lZ + ^inY + Q (n=l,2) 
denote (F^"-parallel lifts of) two complementary nets of /; then the planes (in the quadric Q) 

e" ■.= spaii{F,Z,F",Q} 
of the orthogonal circles c" coincide at corresponding points: 

ef = span{F„ Z,, FJ\Q} = span{F„ Z,, Y„Q}. 

On the other hand, if 

F" = Pijir,) = filZ + fir.Y + Q (n = l,2,3) 

provide three complementary nets of /, then 

Q ec^■.= span{i^" | n = 1, 2, 3} 

for all i e A/, that is, the circles through corresponding points of the three complementary nets 
are, in fact, straight lines in the quadric of constant curvature given by Q. 

Reversing these observations yields four constructions for linear or quadratic conserved quantities 
from geometric configurations of Darboux transforms: 



32) In the smooth case, these orthogonal eheles form a cyclic system, that is, they have a 1-paramctcr family of orthogonal 
sm'faees so that any two orthogonal surfaces are Ribaucour transforms of each other. 
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Thm. 4.11. Let /", n — 1,2, be two Darboux transforms with different parameters fin of a 
discrete isothermic net f so that the circles 

cf^- := span{i^„ F„ /'f } = spa.n{F,,F,,Fp} 

on the edges of M do not coincide for n — 1, 2. Suppose that the /" are antipodal in a suitable 
non-Euclidean space form. Then f has a normalized linear conserved quantity. 

Thm. 4.12. Let /", n ~ 1,2, 3, be three Darboux transforms with different parameters jjLn of a 
discrete isothermic net f so that the circles ' 

cf^- := span{F,,F,,F^"} = spa.n{F„ F, , F^} 

on the edges of M are not cospherical for n = 1, 2, 3. Suppose that the circles 

Ci := spanji;" I n = 1,2,3} 

are straight lines in a suitable space form geometry. Then f has a normalized quadratic conserved 
quantity. 

Proof. We prove both theorems. Thus let A^ = 2 or iV = 3 and let F"-, n = I, . . . ,N denote F'^"- 
parallel light cone lifts of the Darboux transforms /" . Note that the non-degeneracy assumption 
for the edge circles cfj ensures that the F" or F" are linearly independent mod/^ © fj. 

Now let Q e M^''^ \ {0} denote the vector defining the space form. Then 

with suitable functions a„ : Af -^ JR. First, we see that the a„ are constant: 

= dQ.,, - Ell dianhPp + {an)jdF^, = E1i d{ar.).,F^ mod /, ® /„ 
so that d{an)ij = since {F^, . . . , F^) are linearly independent mod fi ® fj. Now, 
= dQ,, = E!Li ^ndF[; = E,li anM„{(F", F), F, - {F", F), F,}, 
showing that 

Z ■■= Ell C^nf^nF" ± / 

defines a sphere congruence enveloped by /: note that |Zp > since the F" and F are linearly 
independent. Hence / has a linear or quadratic conserved quantity if A'^ = 2 or iV = 3, respectively. 



by Lemma 4.9. < 



Note that there are no corresponding theorems for the existence of higher degree polynomial 
conserved quantities for codimension 1 isothermic nets: the non-degeneracy assumption on the 
F"'s, Fi and Fj being linearly independent in K^'^ restricts N to numbers not greater than 3. 

Thm. 4.13. Let /", n = 1,2, be two Darboux transforms with different parameters fin of a 
discrete isothermic net f and let 

M^i^cf :=span{F„Z„i^"} 

denote the orthogonal circle congruences of the Ribaucour pairs {f,f^^), n — 1, 2, with respect to 
an enveloped sphere congruence Z and let 

{ij) K^ Cy := spanjFi, Zi, Fj} = spanjFj, Zj,Fj} 



33) Note that span{Fi, Fj , F"} = span{Fi,Fj. F"} since corresponding edges of a Darboux (or, more generally, Ribaucour 
pair) have concircular endpoints. 

34) That is: Z defines the normal direction at each vertex of the net so that the notion of an "orthogonal circle" is weU 
defined — thus, only the contact element Zi+ RFi at each point i is needed. 
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denote the orthogonal edge circles '. Further, let Q e iR"*'^ dehne a quadric of constant curvature 
so that neither the circles Cij nor the circles ci are straight lines and let 

Cij := spanji^j, Zi,F.j,Q} = spanji^^, Zj,Fj,Q} and ef := spanji^j, Zi, F,, Q} 

denote the corresponding circle planes. Assume that Cij ^ e" for every edge and suppose that the 
circle planes e^ = e^. Then f has a quadratic conserved quantity. 

Thm. 4.14. Let /", n = 1,2, be two Darboux transforms with different parameters /i„ of a 
discrete isothermic net f and let 

i 1-^ cf and {ij) i-^ Cij 



denote the orthogonal circle congruences as in the previous theorem . Assume that Cij ^ c" for 
every edge and suppose that c^ = c^. Then f has a normalized linear conserved quantity. 



The first of these two theore ms is a discrete version of a famous classical theorem, cf. [13 , [14, §84] 



and |[21, Thms. 2.1 fc 2.33] 



Proof. Again, we prove both theorems as the proofs are very similar. 

Let F" denote F^" -parallel lifts of the two Darboux transforms /" and note that Q is F^-parallel. 

First consider the second theorem: as c^ = c^ we have 

F2 = aiF^ +aF + (3Z 

with suitable functions Q;i,a,/?, where ai 7^ since /^ ^ f. Hence 

d{ai),jFl = dPfj -d{aF + (3 Z)^ - {ai)jdFl^ = mod c,,, 

so that Qi is constant since c\ ^ Cij. Moreover, as F^, F^ and F are linearly independent, 

F^ -aiF^ ^aF + (3Z 



defines an enveloped sphere congruence and the claim follows from Lemma 4.9 



In the situation of the first theorem we have 

F^ ^ aoQ + aiF^ +aF + f3Z 
and deduce 

d{ai),,Fl = di^2 -d{aF + p Z\, - d{ao),jQ - {ai)jdFl^ ^ mod e,, , 
so that, again, ai is constant since e\ ^ Cij. Then 

d{ao)ijQ = dP^j - d{a F + /3 Z)ij - (ai)jdFlj = mod Qj, 
showing that ao is constant as well, because Cy was assumed to be a proper circle. Now 

F^ -aiF^ ~ aoQ = a F + /3 Z, 



so that, again, the claim follows from Lemma 4.6. < 

Note that, in the first theorem, we obtain a normalized quadratic conserved quantity as soon as 
we assume that the circles span{F^, F^, F} through corresponding points of /, /^ and /^ do not 
become straight lines in the quadric given by Q. 



35) Note that spanjFi, Z^, Fj} = spaii{Fi, Zj,Fj} by the enveloping condition Dcf. 3.10t 
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5. Discrete cmc nets in space forms 

We are now prepared to define discrete cmc nets in space forms. A smooth isotliermic surface 
/ in the conformal 3-sphere has constant mean curvature _ff in a quadric of constant curvature 
K — — IQp, given by Q € iR^'^, if and only if it has a linear conserved quantity 

P{X) = XZ + Q, 



where Z is the mean curvature sphere congruence ' of /, see [11] or [21, Thm. 2.27]: for smooth 



36) 



isothermic surfaces in the conformal 3-sphere (of codimension 1) it turns out that the top coefficient 
of a polynomial conserved quantity is necessarily its conformal Gauss map ' . Note that, in the 
case of a spherical surface, the conformal Gauss map Z of / is cons tant and / has a constant 
conserved quantity Q = Z (as in the discrete case: see Thm. 3.14 ), hence a linear conserved 
quantity, e.g., 

P{X)=XZ + Q={X + 1)Z. 

In the discrete setting we use this characterization as a definition: 

Def. 5.1. A discrete isothermic net f will be called a discrete cmc net if it is special isothermic 
of type N < I. In particular, if f is special isothermic of type 1 with normalized linear conserved 
quantity 

P{X) = XZ + Q, 

we say that Z is the mean curvature sphere congruence of f in the quadric 

Q = {YeL^\{Y,Q) = -l} 
of constant curvature k — — |QP and that f has (constant) mean curvature 

H:=-{Z,Q). 

Thus discrete cmc nets in space forms are special isothermic net s and a transformation theory is 
readily available to us, cf. [17]: the Backlund transformation, see Def. 4.4 , yields a transformation 
for discrete cmc nets in a given space form preserving the mean curvature, 

Q^Q and H = -{Z,Q) = -{Z,Q) = H 

by Thm. 4.5| , and satisfying Bianchi permutability by Thm. 4.7| ; the Calapso transformation 
provides a Lawson correspondence / i— > /'' for discrete cmc nets, where both the mean and 
ambient curvature change, 



2/ii/ - fi-^ and i/^ 



-(z^Q^ 



H 



M 



since Z^ = T^'Z and Q^ = T^(^Z + Q) by [Thm. 3.13| , but 
remains invariant. 



Def. 5.2. 

dence. 



The Calapso transformation for discrete cmc nets will also be called Lawson correspon- 



The main goal of this section will be to show that our definition generalizes and truly extends 
previous definitions ' from [3|, 16 , [4] and 17 



36) The mean curvature sphere congruence of /, consisting of spheres touching / that have the same mean curvature as 
the surface at the touching points, can be defined using any ambient space form geometry: it can be characterized as the 
conformal Gauss map of /, i.e., the unique enveloped sphere congruence that induces the same conformal structure as /, 
or as the central sphere congruence, i.e., the congruence of spheres that exchange the curvature spheres (via inversion) or, 
equivalently, that have second order contact with the surface in orthogonal directions. 

37) Hence it follows directly that an isothermic surface with linear conserved quantity has constant mean curvature 
H = —{Z, Q) in the space form given by Q. 

38) A change a ^ a = ^ of the cross ratio factorizing funct ion results in a change of the equation |(3.1)| for a linear 
conserved quantity, hence of the hnear conserved quantity, see Lemma 3.6; in particular, Z = Z and Q = cQ. Hence we 
obtain the effect of an ambient homot hety : k = cPk and H = cH. 

39) Note that the mean curvature in |[22]| is defined on the faces of a principal net and therefore different from the mean 
curvature defined here, living on the vertices. 
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5.1. Uniqueness and existence questions 

However, before addressing the relation of our definition with previous approaches we shall discuss 
the construction and uniqueness of linear conserved quantities for a given discrete isothermic net. 

Clearly, given the value of a linear conserved quantity at one point of a discrete cmc net / with 
cross ratio factorizing function a, the linear conserved quantity is uniquely determined as it is a 
parallel section of the isothermic family of (flat) co nnect ions: given an initial value of a linear 
(or polynomial) conserved quantity P{X) we can use (3.1) to determine the values of P{X) at any 
point. Thus, starting with a linear quantity XZ + Q a,t some point of a discrete isothermic net, 
there is a unique parallel section of the isothermic family of connections , but the parallel section 
may fail to be linear or even to be polynomial at other points of the isothermic net — hence we 
may not obtain a linear conserved quantity if we use the "wrong" initial value or if the isothermic 
net is not cmc. 



Using instead the equivalent condition |(3.2)| from [Lemma 3.7| (see also |Cor. 3.8| and |Cor. 3.9 (i-ii)|) , 
we learn that P{X) = XZ + Q is a linear conserved quantity for an isothermic net / with cross 
ratio factorizing function a if and only if 



Q = const, Z J- F and Zj = Zi 



;{{Q,F,)F,^{Q,F,)F,} 



(5.1) 



for any edge (ij) of the domain graph M. Again, we obtain a propagation formula that fixes 
the linear conserved quantity uniquely once an initial value of P{X) is given at some point of a 
discrete isothermic net; however, now the existence of Z as well as the incidence relation Z J- F 
are only satisfied if the isothermic net was, in fact, cmc and if the initial value of P{X) was chosen 
correctly. 

First we address the integrability of the difference equation for the mean curvature sphere con- 
gruence Z, 

Z, - Z, + j^HQ, F,) F, - (Q, F,) F,}. (5.2) 

This equation does clearly not depend on the choice of lift F of the isothermic net /; to simplify 
the computation we may choose a Moutard lift, satisfying |(2.5) , so that 



Zk = 



z^ 
z, 



{Q, Fk - F,) Fj - (Q, Fj) {Fk - F,) 
i^{{Q,F,)Fi~{Q,Fi)F,} 



(5.3) 



on an elementary quadrilateral (ijkl), where we used the Moutard equation (2.6). The last 



expression is symmetric in j and I; hence the propagation equation (5.2) is integrable 
Next we address the incidence relation Z ± F. First consider an edge {ij): from (5.2) 



{Z„F,)^{Z,,F,)+a,,{Q,F, 
Consequently, incidence determines Zi at the center i = i(o,o) 



of a non-spherical vertex star 



if 



we let i(^m,n), rn,n ^ {~li 0, 1}, denote the vertices of a 3 x 3-grid then the equations 



(^., 



.F,„ 



= and (Z. 



*(0,0) 



F,., 



^«(0,0)«(t, 



AQ,F^i„ 



(5.4) 



where rn^ 



1, have a unique solution Zi.^^. since the vertices Fi.^^. and Fi. , of a non- 



spherical vertex star form a basis of K'*'^. Thus an appropriate choice of the initial value Zi.^ ^ 
for the mean curvature sphere congruence at the center of a vertex star ensures that incidence is 
satisfied on all five vertices of the vertex star when using (5.2) to define Z on the corresponding 
3 X 3-grid. At the diagonal vertices Fi , m,n — ±1, we obtain incidence without further 
conditions: let {ijkl) denote an elementary quadrilateral; then, using (5.3)| , (2.6) and [(2.5) , we get 

(Zfc,Ffc) = {Z, + ^^{{Q,F,)Fi-{Q,Fi)F,},F, + ^^{F,-Fi}) 

= {Z„Fi) + ^^-^{(^^ + a,,Q, Fj) - {Z, + a,iQ, Fi)} 

= {Z.,,F,) + ±d^{{Z,,F,) - {Zi,Fi)} 

= 0, 



that is, incidence of Zk and fk- Thus we have proved the following, cf. Cor. 3.4: 



40) Here we use the flatness of the connections in tlie family. 

41) Here we use the fact that any vertex has four neighbours: in a more general quad-graph this argument only works at 
vertices of degree 4, i.e., not at "umbilics" of a discrete net. 
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Lemma 5.3. Let f : {{m,n) \m,n E { — 1,0,1}} -^ S'^ be a non-spherical discrete isothermic 
3 X 3-net and let Q E M^-^ \ {0}. Then f has a unique linear conserved quantity P{X) = XZ + Q. 

Note that, if the vertex star of /(o,o) is cospherical, then the corresponding 3 x 3-net is necessarily 
also cospherical since / is a discrete principal net, i.e., the vertices of its faces are concircular. 



Hence, assuming that / is non-spherical in Lemma 5.3, we have that the vertex star used to define 



Z at its center is non-spherical. Moreover, since any two adjacent vertex stars (i.e., vertex stars 
at the endpoints of an edge) of a discrete isothermic net have two face circles in common, they 
lie on the same sphere if they are cospherical — thus, an isothermic net is either spherical, hence 
type 0, or it has a non-spherical vertex star. 

As a (degenerate) example consider the Moutard lift F of an isothermic 3 x 3-net, i.e., let F satisfy 



{Fi ,Fj) — ttij as in (2.5) , and choose Q so that 



Q^Fn 



(1,0) 



-P'(-l,0)7-P'(0,l) 



-^"(0,-1) I -^"(0,1) 



^(1 



(1.0) 



^ 



{Q,F, 



(7n,n)/ 



Co 



(5.5) 



for TO 



1 and some cq G iR; further let ci :— {Q, i^(o,o)) a-nd observe that, from (2.6) 



{i'{m,n) -^{Ofi),Q) = (f-(„.o).F(o,„)) ^^i"^fi) -^{0,n),Q) 

for m,n = ±1, so that 

{F{m,n),Q) ■ 







Co if m'^ + n? 



C\ if TO 



Now, K5.4)| yields Z(o,o) + coi^(o,o) -L -F(o.o), ■P'(7Ti,n) for vr? + n'^ = 1, that is, Z(o,o) + cqFjo 



orthogonal to a basis of JR ' , hence vanishes. Consequently, (5.2) yields 



0) 



IS 



■'{rrt^n) 



'CQF(^„i,n) if ™^ 



n^ = 0, 



Thus Z \\ F becomes isotropic and yields another Moutard lift of /. Note that, if we choose Q 
so that Co — 0, i.e., so that Q represents the (unique) sphere containing the outer points of the 
vertex star, then Z vanishes for m^ — n^ = 0, that is, we are in the situation that we excluded 



from consideration in Lemma 3.11 , This last situation is a worst-case scenario for a non-spherical 
isothermic net and that, by (5.2) , Z does not vanish as soon as {F, Q) does not, that is, as soon 
as / does not hit the infinity boundary of the space form given by Q. 



Now recall from Lemma 3.11 that, if Z is null but docs not vanish, then it is necessarily a Moutard 
lift of /. In particular, 



Z, 



{711 /n) 



-CoF(,n^n) if m'^ 



-ciF(^,n,7i) if TO^ + n^ = 1 

for some cq, ci S JR and a Moutard lift F satisfying (2.5)| as any two Moutard lifts are related in 
this way. Hence, for to^ + 71^ = 1, 

{F(m,n},Q) ^ --^{Z{rn,n),Q) = COnst, 



SO that we are back in the situation of |(5.5) . Thus avoiding (5.5)| we can ensure that Z becomes 
spacelike and hence a suitable rescaling of the constructed linear conserve d quantity will leave us 



with a normalized linear conserved quantity; hence, according to Dcf. 5.1 



Lemma 5.4. Let f : {{m,n) \m,n G {—1,0,1}} -^ S"^ be a non-spherical discrete isothermic 
3 X 3-Jiet and choose Q S R'^'^ to satisfy, for m,n E { — 1, 0, 1}, 



Q/F(, 



^) 



and g^{F(i 



0) 



F{-i^o),F(^Qi^ 



F{Q,-i),F(^Qi-) 



■^(1,0) I""": 



(5.6) 



wiiere F is a Moutard lift of f. Then f is a discrete cmc net in a space form defined by a suitable 
rescaling of Q. 



Thus, as (5.6)1 imposes only open conditions on Q, any given isothermic 3 x 3-net is cmc in a 
4-parameter family of possible space forms. Enlarging the net and using (5.2) to propagate Z 
will add more incidence conditions, hence conditions on Q; it is therefore natural to expect that a 
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large enough generic isothermic net, will have a unique linear conserved quantity and being cmc 
becomes a condition. In particular, the conditions |(5.1) for a linear conserved quantity on an 
isothermic 5 x 5-net 

/ : {(m,n) I m,rie {-2,-1,0,1,2}}-^ 53 

can be reduced to the incidence conditions on an extended vertex star 

{/(0,0): /(±l,0)i /(0,±l)i/(±2,0),/(0,±2)} 

by the arguments that proved pj^emma 573 : this yields nine linear equations for the linear conserved 
quantity at the center vertex /(o.o) of the net — hence we get existence of a linear conserved 
quantity on any isothermic 5 x 5-net, and we expect uniqueness up to scaling generically, i.e., 
using the scaling freedom to normalize the obtained linear conserved quantity we expect a generic 
isothermic 5 x 5-net to be a discrete cmc net in a unique way. 

However, the following example shows that even "arbitrarily large" isothermic nets can be cmc 
in different space forms — and even with their respective mean curvature sphere congruences 
defining the same contact elements at each vertex. 



For this purpose we recon sider our example (3.5). However, instead of thinking of / as part of a 
"zigzag-plane" as in (3.5) , we now think of it as part of a discrete circular cylinder by letting 

a :^ cos If and /3:=sin(^ 

for some if £ (0, -1), so that ' 

{rjm, °^'^ + ( — 1)" ~'^°^'^ ,n sin ip) = {rjm, cos mp, sin nip). (5-7) 



2 
Further we introduce 



/, 



(m.n) 



{r/m, — cosnip, - 



lip). 



Now let 



Z:=^F*-Q and Q := (1, 0,0, 0, -1) 



(5.8) 



where i^* = ( ^+lfl ,/ 



i-|/-|^ 



denotes again the Euclidean lift of /* , and note that 

(i^,F*)--i|/-r|2^-2, 

so that {XZ + Q,F) = —1. Now observe that |/*P = |/P do not depend on n and hence 



dF, 



(m.n){m-\-l^n) 



d-r{m,n)(m+l,n) — "-^(m,n)(m,ri+l) 



" dF{m,n)(m,n+l) — 0; 



consequently, Z and Q define a linear conserved quantity by Lemma 3.7| for all edges {ij) 



d{\Z + Q), 



Aoi 



{F^.FJ) 



{(XZ + Q, F)jF, - (XZ + Q, F),Fj} = 0, 



where a^ = zt^{Fi, Fj) as in the example (3.5). Since |Zp = 1 the linear conserved quantity is 



normalized, characterizing / as a discrete net of constant mean curvature 

H = -{Z,Q) = i 

in Euclidean space (as jQp = 0) according to Def. 5.1 , as one would expect. 

Superposition of the linear conserved quantities given by (3.6) and by (5.8)| then yields a 1- 
parameter family of normalized linear conserved quantities for /, given by 



Uf* 



-l)^F) - Qo and Q„ := Qo 



I l+cosy n 1 n 1-t-cosy \ 
'yK 2 5 ^7 ^5 ^1 2 ' "^ 



where Qo = (1, 0, 0, 0, —1) and a E M. Thus the isothermic net has constant mean curvatures 

H^ = -(^o,Qo) = i(l + a^) - ^^a 
in the hyperbolic spaces 

of curvatures k„ = — /-, *" — vi- Note that all Z„ define the same contact element at a vertex of 
/, that is, they all define the same "normal direction" at a vertex. 



Finally observe that, as for the discrete net (3.5) with degenerate linear conserved quantity (3.6) , 
the restriction of the domain to { — 1,0,1}^ is again not necessary: t he cir cular cylinder in (5.7) 



can be defined on all of Z,'^ with linear conserved quantity given by |(5.8) . As the two nets, the 



"zigzag-plane" and the circular cylinder, coincide on ^ x {—1,0,1} we obtain an example of an 
isothermic N x 3-net, N G N arbitrary, that is cmc in different space forms. 



42) Remember that m,ne {-1, 0, 1}. 

43) This is the (parallel) Christoffcl transform of /, which is the net of centers of the mean curvature sphere congruence 
in Euclidean space. 
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5.2. Relation with previous approaches 

At this point we are prepared to link the present approacli to discrete cmc nets to previous 
approaches; in particular, we shall discuss how our definition relates to: 



the notion of discrete minimal surface in Euclidean space introduced in [3, Def. 7] 



the notion of discrete cmc net in Euclidean space introduced in [16, Sect. 5] 



the notion of discrete horospherical net (cmc 1 net) in hyperbolic space in [17, Def. 4.3] 



• the definition of a net of constant mean curvature i/ in a space form of curvature k, where 
H^ + K>Q, suggested in ^Tt} . 

In this context we shall also discuss the relation of our notion of "mean curvature sphere" with 
that of |[4, Sect. 4.5]| and the notion of "central sphere congruence" in [7, Sect. 3]. 



First we wish to make contact with the definition of a discrete cmc net in Euclidean space given 
[16] : recall that a discrete isothermic net / : M^ -^ R^ is called cmc in [16] if it has a parallel 



(isothermic) net, i.e., a simultaneous Christoffel and Darboux transform /*. The (constant) 
distance of this parallel net yields (up to sign) the mean curvature of both nets, 

TT 1 

-^ \f-f\- 



Note that any constant multiple of {df, df*) is a cross ratio factorizing function , so that, without 
loss of generality. 



fly = -f {df^j,df*j), that is, d/*- = 
Further, since [/* — /[^ = const, (2.2)| yields 

(r-/).,±dA„c^, 



H\df,j\- 



■df,. 



that is, iJ (/* — /) defines a (unit) normal field for / as well as for /* in the sense of [22] , cf. [16] 
Consequently, 



d{\m,=2{n,dn)^ 



H\df, 



-\^{fij,dfij) - -TJ\dj^d{\f\ )ij. 



Now suppose / : IvP -^ ]R^ is a discrete cmc net in this sense, with parallel cmc net /* and 
(constant) mean curvature H. Let F and F* denote the respective Euclidean lifts and let 

Z:=HF*-^Q and Q == (1,0, 0,0, -1) 



as in the above example of a discrete circular cylinder, see (5. 

{z,F)^-§\r~f\^ + ^ 

and 



Then 



2H 



dZi. 



°'3 rJP.- 



HdF* 






0, 



so that \Z + Q defines a (normalized) linear conserved quantity of / by Lemma 3.7. Moreover, 

K = -[Q[2^0 and -{Z,Q)=H. 

Conversely, suppose that / is a discrete isothermic net with cross ratio factorizing function a and 
normalized linear conserved quantity XZ + Q so that 

K = -|<9I^ = and H = -{Z,Q) ^Q. 



44) Given a cross ratio factorizing function a of an isotlicrmic net one defines the Christoffel transform /* by 

-1 _ _JHj ir. .. 



Ci ='^ii{dft3) 



M/iiP 



as a is only defined up to constant multiples, so is df* . In the case of the parallel net of a discrete cmc net, however, there 
is a canonical scaling for /* . 
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F* := ±{Z+^Q) = ^{2HZ + Q) 



denote the complementary net of /, see Def. 4.8 . Clearly, F* defines a Darboux transform of /, 

T^^F* = const. 



and {F\Q) = -1. Choosing a Euclidean lift F of /, i.e., (F, Q) = -1, |(3.2)| yields 



dF*, = 



■ dF,,, 



ij ~ H{F,,Fj) ""^ y 

Now let, without loss of generality, Q = (1, 0, 0, 0, —1) so that 

F* = (i±^,r,i^) and F=C-±\l^,f,'-]I^). 



(5.9) 



^^^^dU 



Then |(5.9)| gives 

Hence /* is also the Christoffel transform of / ; NP -^ ]R^, so that / is a cmc net with parallel 
cmc net /* in the sense of |[16]| . 

Thus we have proved the following 



Thm. 5.5. A discrete isothermic net in M^ is cmc in Euclidean space in the sense of [16] if and 
only if it is cmc in the sense of Def 5. 1 . 



In [4, Sect. 4.5] , a notion of mean curvature sphere for a discrete isothermic net in Euclidean space 
is introduced. We shall see that this mean curvature sphere is the same as our mean curvature 
sphere Z in the case of a discrete cmc net. 

First recall [[4, Def. 12] : given a vertex star {fi,^ „, € R^ \ m^ + n^ < 1} of a discrete isothermic 
net with constant (negative) cross ratio function qijki = -^ there is a unique point Ci,,-, p. so that 



\ji(i,o) '^i(o,o)\ ~ \ji(-i,o) ''*(o,o) I ^^" \ji(o.i) '^i(o,o)\ ~ \j^i(o.-i) ''«(o,o) I (5.10) 



and 



l-''(l,0) '^'(0,0)1 l-''(0.0) "^'(0,0) I l-''(0,l) "^'(0,0)1 l-''(0,0) ^'(0.0)1 



°'(0,0)»(1,0) "'(0,0)'(0,1) 

This point Ci.g „, is the center of the mean curvature sphere at fi.g g, and its radius 



Ho.o) 



\ji{o,o) '^i(o,o)\- 



(5.11) 



(5.12) 



Now suppose that / is a discrete cmc n et in Euclidean space, with Q as above and mean curvature 
sphere congruence Z, and go back to (5.4) : let F denote the Euclidean lift for / as before and 
write 

Z=i(l + |c|2~r2,c,l-|c|2+r2) 

in terms of its center c and radius r; then 

0= \^i(p,a)T^i{0,0)/ = ^2r,^„^j (1/^(0,0) ~ ^i(a,o)\ ~ ^i(o,o) ^ 



is equivalent to (5.12), while the remaining four equations of (5.4) read 



<^-(o.o)'^'(....)) 



l■''(m,,^) '•'(0,0)1 ''(0,( 



'(0,0)»(m,»i) 



clearly implying (5.11) ; the equations (5.10) follow since the cross ratio function is constant, so 
that 

'^«(o,o)«(i,o) ~ '^■'(o,o)«(-i,a) ^^"- '^«(o,o)«(o,i) ~ '^■'(o,o)Ho,-i) ■ 

Hence we have proved: 
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Thm. 5.6. The mean curvature sphere Z of a discrete cmc net f : M -^ ]R^ hi Euchdean space 



in the sense of Dcf 5.1 is the mean curvature sphere of f in the sense of [4, Sect. 4.5 



In fact, we have seen slightly more: the equations (5.4) define the mean curvature sphere of |4| 
for any isothermic net in Euclidean space with constant cross ratio function ' . 

This suggests to use 1(5.4) to define the mean curvature sphere of an isothermic net in any space 
form. Note that these mean curvature spheres of an isothermic net are, in contrast to the smooth 
case, not Mobius invariant: different choices of the ambient space form given by Q will, in general, 
lead to different mean curvature spheres Z . In particular, this also shows that the mean curvature 
sphere defined in this way is generally different from the central sphere of [7] , which is Mobius 



invariantly related to the isothermic net as it is characterized by incidence only. 



Following ideas from [17], the Lawson correspondence of Def. 5.2 can be used to define discrete 
cmc nets in space forms with 

as Calapso transforms of cmc nets in Euclidean space. Using permutability theorems the above 
characterization of cmc nets in Euclidean space by the existence of a simultaneous Darboux and 
Christoffel transform, i.e., the existence of a parallel cmc net can then be carried over to other 
space forms to obtain an alternative characterization, similarly to the way in which horospherical 
nets in hyperbolic space can be characterized in two ways, see [17, Lemma 4.2] or [18, §5.7.37] 



Namely, let /* denote the parallel cmc net of a discrete cmc net / in Euclidean space and consider 
their Calapso transforms /^ and (f*)^ — which are only determined up to Mobius transformation. 
Then /'*' and (/*)'*' can be positioned in S^ so that: 

(i) they form a D ar boux pair wi th parameter —A since / and /* form a Christoffel pair, see 



[17, Cor. 3.24]| or |[18, §5.7.34] 



(ii) they form a Darboux pair with parameter /i — A since / and /* form a Darboux pair with 
some parameter ^, see [17, Cor. 3.27] or [18, §5.7.35| . 



That is, there are two ways to position (f*)^ in S^ so that it is a Darboux transform of /'^ with 
different parameters or, otherwise said, f^ has a pair of Mobius equivalent Darboux transforms. 



More precisely, given a Calapso transform of a discrete cmc net in Euclidean space, (4.13) provides 
two antipodal Darboux transforms: the antipodal map identifies the ambient space form. 

provides the precise formulation of the converse: normalizing the linear conserved 



Thm. 4.11 



quantity from the proof. 



(1 + A^i)aiFi + (1 + A^2)a2^^ 



we obtain |P(A)|^ = 



2 _ (l+AAti)(l + AM2) 



^1/^2 



as its squared norm, hence 



H^ 



V 2fJ,ifJ,2 ' MlM2 






since \x\ 7^/12, so that the isothermic net / is indeed a Calapso transform of a cmc net in Euclidean 
space. 

Finally we discuss minimal nets in Euclidean space and horospherical nets in hyperbolic space, 
that is, the discrete constant mean curvature nets with 



H' 



K = 0. 



45) Besides the sphere, the equations |^5.4)| also determine the scaling of its representative in Minkowski space. 



'22 



46) Recently, a new approach has come into focus, where the "mean curvature" of a discrete p 
as a function on faces, via the area change of a face when varying through parallel nets, see 
leads to the same class of discrete minimal or constant mean curvature nets as the one discussed here 

47) Using aMoutard lift satisfying (2.5), as for (5.5), ^2.6) and (5.4) yield 



inci pa] net in JR^ is defined, 



and [S]. This approach 
9^ 



'■^'(0,0) '-'^'(t, 



,) = o 



'''(0,0)'(m,0) 'Q'-'^Vm.O)) ~ "icO.O) '(0,n) \Q^^i(0,n)/ 



for m,n = ±1. Hence the "mean cmvature sphere" defined by (5.4) is the "central sphere" of |7| (cf Lemma 2.2) if and 
only if, for vn? -\- v? = 1, 



{Q,ai^^ 



0,0)\rr,,-n)^\r, 



.,) = 



const. 
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First recall that horosphcrical nets in hyperbolic space can be defined as Darboux transforms of 
their hyperbolic Gauss maps [17, Def. 4.3], that is, as Darboux transforms of a spherical net: 
hence Lemma 4.10| identifies them as the discrete cmc nets with 



H' 



K = and K 7^ 



in the sense of Def. 5.1. On the other hand, horosphcrical nets can equivalently be characterized as 



Calapso transforms of a discrete minimal net in Euclidean space in the sense of [3] and, conversely, 
any discrete minimal net in Euclidean space gives rise to a Lawson family of horospherical nets, 
17, Lemma 4.2 . Hence we can reverse the argument to conclude that the discrete minimal 



nets in the sense of [3| are those discrete cmc nets with 



H' 



K = and K = 



in the sense of Def. 5.1 . Note that [3, Thm. 8] provides two equivalent characterizations of discrete 
minimal nets in Euclidean space: 

(i) as Christoffel transforms of spherical isothermic nets (their "Gauss maps" ) and 

(ii) by the fact that their mean curvature sphere congruence consists of planes . 

We summarize these discussions in the following 

Thm. 5.7. A discrete isothermic net is 
(i) minimal in ]R^ in the sense of [3j iff it is cmc with H =^ 



in the sense of Def. 5.1 ; 



(ii) horospherical in the sense of ¥171 iff it is cmc with —H = k < in the sense of Def. 5.1 



5.3. Discrete cmc surfaces of revolution 

We conclude by discussing our example of a discrete surface of revolution in more detail: here we 
will be interested in the construction of discrete cmc surfaces of revolution with prescribed mean 
curvature in a given ambient space form. 

Thus we consider a discrete surface of revolution 



(^m,n) ^ F(„,„) = (-l)™(Af„ + $„C) £ JR^.i ^ ^^^2 



see 



2.8), with cross ratio factorizing function 



Uij = a{F^,F.j 



(5.13) 



(5.14) 



since F is a Moutard lift . By Cor. 3.1(j a linear conserved quantity P{\) = XZ + Q has the same 
rotational symmetry as the net if and only if the vector Q defining the ambient space form has this 
symmetry; in particular, the linear conserved quantity of an equivariant cmc net is rotationally 
symmetric with 

P(m,n)W - P,ilW + (-l)™+iAp,„(A)$„C = X{Zi - a(Q, A/™) <D„C} + Q, 

(3.8)1 and Lemma 3.6. Hence 



where Z-^,Q e JR ■ , see 



^(?Ti,n) 



a{Q,M„ 



Oi{QiF{m,n+l)) F{m,n), 



(5.15) 



48) Assuming, as in [3] that a.-, = ±5 a very similar computation as the one leading to Thm. S.q shows that our statement 
here is equivalent to [3. Def. 7] : with the Euclidean lift F and Z = (d. n. —d) . where n is the unit normal of the plane 
described by Z and d = («, /) its distance from the origin, the equations (5.4) become 

\"»(0,0) ' f'{m 



^'(0,0)/ ~ ^^'{0.0) '^Hm.,!)/ - °'\0,0)H 



= ±5. 



49) We shall use the scaling freed nnn n iii th 



quantity that we will construct, see Lemma 3.6 



cross ratio factorizing function later to normalize the linear conserved 
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where S,,,, — S, 



(rn,n)(rn,n+l) 



is the curvature sphere (3.4) on an edge in the rotational direction: 



note that this family of curvature spheres does only depend on m, so that our discrete surface of 



revolution is the envelope of a 1-parameter family of spheres in the sense of Def. 3.10 



Our first aim is to formulate the condition (3.2) for a linear conserved quantity of the given form: 
clearly 

0=(^,F)-(-l)"(^,M), 

that is, the incidence relation is again expressed by orthogonality; further, 

= rf-Z'(,n,n)(m+l,n) + Ct{{Q, F{m,n)) F(m + l,n} ~ {Q , F{m+l,n)) F(m,n}} 
= dSrn,m+l ^ 2a{Q, Mra.m+l) dMm,rn+l 

and 

— dZ(rn,n)(m,n+l) + <^{{Q , F{m,n)) F(m,n+1) — {Q : F{m,n+1)) -F'(m,n)} 

are identically satisfied. Hence we obtain: 

Lemma 5.8. Let Q € JR^'^ and 2'(^„) := Sm — ct{Q,F(^m,n))F(^m,n)> where F is a Moutard lift 
(5.13) of a discrete surface of revolution and m i— > Sm S JR^'^ is a discrete 1-parameter family of 
spheres. Then P{X) := XZ + Q defines a linear conserved quantity for f with respect to (5.14, as 
a cross ratio factorizing function if and only if 



= {S, M) and dS = 2a(Q, M) dM. 



(5.16) 



Note that the equations (5.16) determine Sm and Sm+i (and, hence, Zm and Zm+i) up to a 
common orthogonal offset from the plane spanned by Mm and Mm+i- Also prescribing the 
mean curvature H , 

(5, g> = {Z, Q) + a{M, Q)2 = -H + a{M, Q)^, (5.17) 



51) 



we obtain six equations' 
{Sm,Q) =-H + a{Mm,Qy 

{Sm,Mm) =0 



{Sm+1,Q) =-H + a{Mm+iM? 

{Sm+l,Mm) = —OL{Q,Mm,m.+ l) \dMm,m+l\ 



{Sm,Mm+l) = -a{Q, Mm,m+l) \dMm,m+l? {S 



m+1, Mm+1 



= 



that determine Sm and Sm+i as soon as (Q, Mm, Mm+i) is a basis of K^'^. 

In order to understand this restriction geometrically we analyze what happens when Mm, Mm+i 
and Q are linearly dependent, that is, when Q is in the (1, l)-plane spanned by Mm and Mm+i . 
First recall that the light cone in ]B?'^ is the axis of our discrete surface of revolution (2.8) ; 
similarly, for fixed n, 

span{i^(„j_„),f(„+i^„),$„C} = span{M,„,Af„+i,$„C} =: Cm,m+i 



defines the circle through Fim.n) and i^( 



{in+l.n) 



that intersects the axis orthogonally — which 



becomes a geodesic in the quadric Q of constant curvature (3.3) given by Q as soon as 



^ t Cm.ni+l- 

Thus prescribing an edge of a meridian curve. Mm and Mm+i, and a space form, Q, so that 
the orthogonal circle Cm.m+i of the axis passing through F(m,n) and F(^m+i,n) does not bec ome a 



straight line, any choice of H will lead to a unique solution Sm and Sm+i of the equations (5.16) 



and (5.17) , hence to a unique linear conserved quantity for the discrete net obtained by rotating 
the edge. 



The to p coefficient Z oi a. linear conserved quantity constructed fr om a solutio n S of (5.16) and 
(5.17) is, with S, spacelike since S ^- M and has constant length by Lemma 3.5 . The idea is then 



to use our scaling freedom a in the cross ratio factorizing function (5.14) to obtain a normalized 



50) Recall that M in [2.8) takes values in one c ompon ent of t he hyp erbolic quadric |Fp = —1 in j R^'-^. 



51) Note that d{{M,Qf) = 2{M,Q){dM,Q) by ;2.2) so that (5.17) 



only adds one real equation to (5.16) 
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linear conserved quantity, as sought in Def. 5.1: note that, in contrast to Lemma 3 .61 , where only 
the conserved quantity condition played a role, the prescribed mean curvature equation (5.17)| 
causes a more complicated dependence of S" on a. In particular , 



\S\-=<^H^-i{a+^Hf 



where we let 

A 
C 
B 
A 



|gAA//™AAf„+i|2, 

\dMrn.m+l\ (Afm,m+1,Q), 



-\dMm,m+l\ 
' |M„AA//,„+iP 



{|M„,™+i|2(Af,„,Q)(M„+i,0>+2(M„,,„,+i,Q)2}, 



Note that \Q A M„ A M^+ip, \Mm A M^+ip < 0. Also, the assumption M^, M^+i / Q that 
the endpoints of our meridian edge do not lie in the infinity boundary of the space form defined 
by Q implies that B and C do not simultaneously vanish; as a consequence ^ > 0. 

Thus the equation \S\'^ = 1 can be solved for a if and only if C^H^ < A. 

Moreover, as a is a factor of our cross ratio factorizing function, we seek a non-zero solution of 
the equation: since A = C^H^ clearly implies H ^ we need to exclude 



A - C^H^ = B = => ij2 = 



Thus the equation \S\'^ — 1 has a non-zero solution a if and only if 
C^H^ < A and C^H"^ = A^H^^ 



A 



|Af,„AM„,+i|2 



(5.18) 



Note that, in the case of strict inequality in |(5.18)| , C^H^ < A, our construction will, in general, 
provide two different linear conserved quantities . However, in the case iJ = of a minimal 
surface, the equation jS'p = 1 has always exactly one non-zero solution a. 

Lemma 5.9. Prescribing an ambient space form and a mean curvature, the equivariant discrete 
surface of revolution obtained by rotating a single edge in tlie space form lias a normalized linear 
conserved quantity as soon as: 

(i) the orthogonal circle of the axis of revolution passing through the endpoints of the edge is 
not a straight line in the ambient constant curvature geometry, and 

(ii) the mean cur vature is not chosen too large; more precisely, the mean curvature H satisfies 
the constraint (5.18, . 



Having equipped an initial edge of a meridian curve for a discrete equivariant cmc net with 
prescribed mean curvature H and ambient space form Q with a suitable mean curvature sphere 
at both endpoints, we shall now investigate how to propagate the meridian curve M and its 
enveloped 1-parametcr family of spheres S to "build" a larger equivariant cmc net. That is, we 
aim to construct Mm+i and Sm+i from the data at the other endpoint of the edge. Mm, Sm, Q- 

As a discrete analogue of a constant speed parametrization in the hyperbolic plane of the meridian 
curve we prescribe a constant cross ratio factorizing function a along the meridian curve so that 



(Af, 



m+lj 



M„ 



(1 



-(1 



for all TO and some c G JR, see (5.14) . Note that, since we seek Mm and Mm+i to take values in 
the same hyperbolic plane, 

12 >0. 



£. — L\f]M 



m,m+l I 



52) Note that the symmetry of the formula confirms that |5mP = \Sm+i\^ ■ 

53) As a consequence, our discrete net will have a polynomia l conserved quantity of degree by |Cor. 3. 



riun. 3.1 



] 



- which is not 



too surprising since the net is clearly spherical, cf 

54) Recall that a hyperbolic constant speed parametrization of the meridian curve of a (smooth) surface of revolution 

leads to a conformal curvature line parametrization of the surface up to constant rescaling of the parameters. 
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A second equation obtained from (5.16), 

= {M„,+i,S,n) + a{Q, M„,„+i)|dM„,„,+i|2 = {M,r,+i,S,n + cQ) + c(Af„, Q), 

then confines Mm+i to a line in JR^'^ as soon as Mm and Sm + cQ are linearly independent or, 
equivalently, as soon as 

X,„ := 5„ + c (Q + (Q, M„) M„) ^ 0. 

Observe that, requiring the meridian curve to not cross the axis of rotation into the other hyper- 
bolic half plane nor to hit the infinity boundary Q^ of the ambient space form \ 

^ c (Q, M„,){(Mm+i, M„> - 1} = (M„+i, X„), 

so that Xjn 7^ describes a sphere intersecting the axis of rotation orthogonally and containing 
the point of the meridian curve given by Mm since Xm -L M„i- In particular, X ^ at both 
endpoints of a "proper" meridian curve edge. 

Thus 

T^,f -SO'+'^M {(a+cf-a^){Q,M,^) X^ ^ . , Y^ 

for a suitable t E R, where Ym is orthogonal to Xm in the (Euclidean) plane Mm in JR^^^ and has 
the same length, 

Ym-LMm,Xm and \Ym\^ = \Xm\^. 



As we wish Af to take values in the hyperbolic plane 



= 1 + lA/, 



m+l| 



ix„r 



K^^ 



^}, 



where we have used (5.17), so that we obtain two candidates for Mm+i as soon as 

l-2cH -c^K>0. (5.19) 

When propagating the meridian curve with constant cross ratio factorizing function, then one of 
the two solutions must give the predecessor Mm-i of Mm', hence the propagation of a meridian 
curve is unique if it is possible. Observe that the condition [(5.19) does not depend on m; hence it 
is automatically satisfied as soon as a "seed" meridian curve contains more than one edge ' 



Once Mm+i is constructed (5.16) yields 

Sm+l '■— Sm + 2a{Q, Mm.^m+l) dMm,m+l', 

it is then straightforward to verify incidence 

{Sm+uMm,+ l) = {Mm+l,Sm + cQ) + c(Q, A/,„) = 



and the mean curvature (5.17) being constant, 

{Q,Sm+i)-a{Q,Mm+i) = {Q,Sm)-a{Q,Mm). 



Hence, by Lemma 5.8, we have succeeded in propagating the meridian curve of a discrete cmc net 
of revolution: 



55) Wc do, however, not exclude the possibihty of the meridian curve crossing the infinity boundary. 

56) In ease we have one edge of a seed meridian curve we infer that the quadratic equation has at least one solution, i.e., 

l-2cH -c^K> 0. 

In the ease 1 — 2cH — c^k = of one solution, the described propagation procediu'c will then alternate the two endpoints 
of one edge. 

57) Recall that the equations '5.16) also ensure that |Sp is constant. 
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Lemma 5.10. Let a point and unit normal of the meridian curve of a discrete cmc net of revo- 
lution in a space form be given. Prescribing a (constant) cross ratio factorizing function c and a 
Moutard factor a, the meridian curve can be propagated uniquely in either of two directions to 
obtain a discrete cmc surface of revolution as long as: 

(i) — > and the point of the meridian curve does not lie in the inhnity boundary of the space 
form, and 



(ii) 1 — 2cH — (?K > 0, where the mean curvature H is given by (5. 1 7 , and k is the ambient 
curvature. 




Fig. 2. Discrete cmc torus in S'^ 



These two lemmas now provide a method of construction for discrete cmc nets of revolution in a 
prescribed ambient space form and with prescribed mean curvature: 

Construction. Choose an ambient space form Q € iR^'^ and a mean curvature H E JR. 

Choose an initial edge Mq, Mi £ H^ C JR^'^ so that: 

(i) {Q,Mo), {Q,Mi) ^ 0, that is, Mq and Mi do not lie in the inhnity boundary of the space 
form ' Q; 

(ii) (Q,Mo,Mi) is a basis of IR"^'^, that is, the straight line in H^ through Mq and Mi is not 
straight in the chosen space form Q; and 



(iii) the constraint (5.18^ is satisfied with the chosen mean and ambient curvatures. 

Next, construct the sphere s So , Si £ S*^ C B?'^ to s atisfy the conserved quantity and mean 
curvature conditions (5.16) and (5.171 , see Lemma 5.£ ; generically, there is a choice of two such 
sets of spheres. 

Then, propagate this "seed" meridian curve using a constant cross ratio factorizing function, see 



Lemma 5.10: this yields a unique and proper propagation of the meridian curve if (5.19^ is satisfied. 



A variant of this construction was used to obtain the discrete cmc torus in S'^ shown in Fig. 2 



as well as to construct the discrete minimal net in two copies of H^ shown in Fig. 1 : the plane 



shown in the figure indicates the common infinity boundary of the two copies of the ambient H^. 



58) Note that this condition is sufficient but not necessary, as Fig. 1 suggests. 
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